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Preface 


The  Air  Force  Weapons  Laboratory  (AFWL/ALO)  expressed 
a need  to  know  what  characteristics  of  reflected  laser 
radiation  could  effect  the  operation  of  the  Adaptive  Laser 
Optical  Techniques  (ALOT)  tracking  and  pointing  system. 
Because  the  radiation  that  is  used  for  tracking  is  reflected 
from  rough  surfaces,  the  far  field  intensity  pattern  has  a 
complex  distribution  of  irradiance.  This  complex  distri- 
bution of  irradiance  is  called  speckle. 

A computer  model  to  calculate  and  plot  laser  speckle 
patterns  was  developed  and  tested.  Because  a rough  surface 
scatters  laser  light  in  all  directicms,  predicting  the  ampli- 
tude and  phase  of  the  electric  field  of  the  scattered  light 
is  very  difficult.  A complex  aperture  function  model,  used 
to  describe  this  electric  field,  was  incorporated  into  the 
speckle  model.  This  complex  aperture  function  model  enables 
the  user  to  input  to  the  program  the  surface  height  measure- 
ments, the  surface  reflection  coefficient  distribution,  and 
the  incident  laser  radiation  profile  at  the  scattering 
surface. 

The  program  was  developed  to  run  on  the  CDC  6600  com- 
puter. It  was  tested  with  simulated  surface  height  data  and 
compared  to  experimental  observations.  The  results  obtained 
from  the  program  compared  very  favorably  with  theory.  The 
results  of  this  thesis  and  further  application  of  the 
speckle  program  should  be  very  valuable  in  understanding  the 


ii 


effects  of  speckle  patterns  on  the  operation  of  the  ALOT 
system. 


I appreciate  the  counseling  of  my  advisor.  Major  Richard 
Potter,  and  the  assistance  of  my  wife,  Karen,  in  the  prepar- 
ation and  typing  of  this  report. 

William  E.  Landis 
Captain,  USAF 
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Abstract 

Speckle  patterns  produced  by  laser  radiation  scattering 
from  rough  surfaces  can  be  calculated  and  plotted  using  the 
SPECKLE  program  developed  in  this  report.  The  scattering 
surface  characteristics  that  are  required  by  the  program 
are;  surface  height  measurements,  surface  reflection  coef- 
ficient distribution,  and  incident  laser  radiation  profile. 
These  three  inputs  are  required  to  calculate  the  complex 
aperture  function  at  the  scattering  surface,  which  is 
Fourier  transformed  to  produce  the  speckle  pattern. 

Through  experimental  investigation  and  program  verifi- 
cation it  was  determined  that  the  typical  speckle  width  in 
a pattern  is  directly  related  to  range  and  laser  wavelength, 
and  inversely  related  to  scattering  surface  cross  section. 

The  highest  spatial  frequency  present  in  a speckle  pattern  is 
proportional  to  the  inverse  of  the  typical  speckle  width. 

The  report  also  showed  that  the  randomness  of  the  scat- 
tering surface  reflection  coefficient  also  effects  the 
speckle  pattern.  This  is  an  important  concept  in  camouflage 
and  very  important  in  the  investigation  of  the  ALOT  system 
performance . 
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MODEL  FOR  PREDICTING  LASER  } 

SPECKLE  PATTERNS  RESULTING 
FROM  ROUGH  SURFACE  SCATTERING 


I.  Introduction 

The  design  of  an  optimum  laser  weapon  system  must 
include  provisions  for  pointing,  focusing,  and  holding  the 
laser  beam  on  target  accurately  enough  to  provide  sufficient 
damage.  Conventional  tracking  and  pointing  systems  using 
radar,  infrared,  or  visual  techniques  do  not  provide  the 
required  accuracy.  A concept  for  fine  pointing  and  auto- 
matic focusing  of  a laser  beam,  referred  to  as  the  Adaptive 
Laser  Optical  Techniques  (ALOT)  system,  is  now  being  investi- 
gated at  the  Air  Force  Weapons  Laboratory  (AWFL) . The  ALOT 
system  measures  the  intensity  of  the  reflected  laser  radia- 
tion from  the  target  and  then  generates  a set  of  corrective 
pointing  control  signals  to  maximize  the  return.  The  ALOT 
system  has  performed  well  throughout  a series  of  laboratory 
experiments.  However,  the  experiments  were  run  under  near 
ideal  atmospheric  conditions,  and  the  effects  of  absorption 
thermal  blooming,  and  turbulence  were  not  rigorously  inves- 
tigated. In  addition,  the  phenomenon  of  speckle,  a result 
of  rough  surface  reflection,  was  observed  but  its  overall 
effect  on  the  ALOT  system  was  not  theoretically  investigated. 
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The  goal  of  this  thesis  is  to  investigate  rough  surface 
scattering  of  a reflected  laser  beam,  excluding  atmospheric 
deteriorations,  and  to  develop  a computer  model  that  calcu- 
lates and  plots  a far  field  speckle  pattern  for  a given 
reflecting  surface. 

Origins  of  Speckle 

Operation  of  the  first  continuous  wave  helium  neon 
(CW  HeNe)  laser  in  1960  revealed  that  objects  viewed  in  a 
highly  coherent  light  acquire  a speckled  appearance.  The 
physical  origin  of  the  observed  speckle  phenomenon  was  quick- 
ly recognized  by  the  early  workers  in  the  laser  field  (Ref.l; 
Ref. 2).  When  monochromatic  light  is  reflected  from  a rough 
surface  the  optical  wavefront  resulting  at  a distant  point  is 
the  summation  of  many  dephased  wavefronts,  each  produced  by 
a different  surface  scattering  element.  The  distances  tra- 
veled by  the  numerous  wavefronts  from  a point  on  the  scat- 
tering surface  to  a unique  point  in  the  observation  plane, 
may  differ  by  several  integer  and  noninteger  wavelengths  of 
the  reflected  radiation.  Interference  of  the  dephased  wave- 
fronts  produce  the  speckle  pattern. 

The  random  interference  phenomenon  of  laser  speckle 
parallels  other  branches  of  physics  and  engineering.  Some 
of  the  more  modern  direct  analogies  are  coherent  imagery  in 
synthetic  aperture  radar  (Ref. 3),  and  radio  wave  propagation 
(Ref .4; 188-267,  Ref . 5) . 
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Theoretical  Investigation 

The  overall  goal  of  this  thesis  was  to  develop  a com- 
puter model  to  calculate  and  plot  speckle  patterns  resulting 
when  laser  energy  is  reflected  from  rough  surfaces.  Before 
the  computer  program  was  developed,  an  extensive  literature 
search  was  made  to  determine  what  theories  have  been  hypo- 
thesized and  developed  concerning  the  prediction  of  speckle 
patterns.  The  bibliography  and  supplemental  bibliography  of 
this  thesis  lists  numerous  books  and  articles  concerned  with 
the  speckle  phenomenon. 

The  initial  paragraphs  of  the  theory  section  outline 
the  history  of  the  investigation  of  speckle.  The  well  known 
Huygens-Fresnel  diffraction  formula  is  the  basic  equation 
used  in  the  development  of  the  computer  model.  Although  the 
Huygens-Fresnel  diffraction  formula  can  be  formulated  from 
experimental  observations,  it  can  also  be  derived  from  the 
classical  scalar  wave  equation. 

From  the  Huygens-Fresnel  equation  it  is  shown  that  the 
far  field  diffraction  pattern  is  the  two  dimensional  Fourier 
transform  of  the  scattering  aperture.  As  a review,  the  last 
part  of  the  theory  section  discusses  discrete  data  process- 
ing, including  sampling  theory  and  the  fast  Fourier  transform. 

Experimental  Investigation 

Before  any  attempt  was  made  to  develop  a computer  model 
that  would  calculate  speckle  patterns,  several  laboratory 
laser  experiments  were  accomplished  to  observe  speckle 
patterns.  The  two  properties  that  were  investigated  were 
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speckle  size  variations  related  to  range  variations  and 
speckle  size  variations  related  to  reflecting  aperture  size 
variations . 

Direct  observations  and  computer  analysis  of  photomul- 
tiplier data  were  used  to  analyze  the  experimental  results. 
The  initial  method  used  to  observe  and  analyze  the  speckle 
patterns  was  to  record  a portion  of  each  speckle  pattern  on 
film.  Photographs  provided  an  easy  means  of  observing 
trends  in  the  speckle  patterns  as  range  and  aperture  size 
were  varied,  but  for  statistical  analysis  the  photographs 
were  a poor  medium  of  recording.  To  provide  a higher  reso- 
lution of  contrast  across  the  speckle  pattern,  a photomul- 
tiplier was  used  to  measure  the  reflected  intensity  at  the 
observation  plane.  The  photomultiplier  data  was  then  ana- 
lyzed, by  a fast  Fourier  transform  computer  program,  to 
obtain  the  spatial  frequency  content  of  the  speckle  pattern. 
This  transformed  data  provided  the  best  method  of  observing 
the  speckle  pattern  characteristics  as  predicted  by  the 
Huygens-Fresnel  diffraction  formula. 

Computer  Model  Used  to  Calculate  Speckle  Patterns 

After  an  extensive  literature  search  on  the  phenomenon 
of  speckle  and  an  experimental  investigation,  the  computer 
model  was  developed.  Basically,  the  computer  model  performs 
the  two  dimensional  Fourier  transform  of  the  complex  ap2r- 
ture  function.  The  complex  aperture  function  describes,  in 
detail,  the  reflected  field  distribution  immediately  in 
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front  of  the  reflecting  surface.  For  a smooth  reflection 
surface,  the  complex  aperture  function  is  a constant  value 
across  the  aperture,  however,  for  rough  surface  reflection 
there  are  amplitude  and  phase  variations  across  the  aperture. 
The  model  used  to  predict  the  complex  aperture  function  was 
suggested  by  Goodman  (Ref. 6)  in  his  study  of  the  statistical 
properties  of  laser  speckle  patterns. 

The  chapter  describing  the  speckle  computer  model  de- 
tails the  formats  used  for  the  computer  data  and  describes 
the  possible  outputs  from  the  program.  The  input  used  to 
calculate  the  complex  aperture  function  is  an  array  of  sur- 
face height  measurements  across  the  scattering  surface.  To 
accurately  describe  a scattering  surface,  an  enormous  number 
of  height  measurements  are  required,  however,  the  computer 
core  available  to  calculate  the  two  dimensional  Fourier 
transform  is  limited.  This  limited  core  size,  which  re- 
stricts the  number  of  input  data  samples,  is  probably  the 
major  deficiency  in  the  computer  model. 

Analysis  of  the  Speckle  Computer  Model 

A ma3or  portion  of  this  thesis  is  devoted  to  the  pro- 
cedure used  in  the  verification  of  the  computer  model  used 
to  calculate  speckle  patterns.  Various  complex  aperture 
shapes,  with  known  two  dimensional  Fourier  transforms,  were 
input  to  the  speckle  program  and  the  results  compared  with 
theory.  Without  exception,  the  computer  model  calculated 
and  plotted  speckle  patterns  that  compared  very  favorably 
with  theory.  When  the  apertures  were  modeled  as  randomly 
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rough  reflecting  surfaces,  the  speckle  patterns  could  not  be  ' 

compared  directly  to  the  known  results,  however,  several 
characteristic  features  were  observed  and  compared  to  theory. 

The  two  major  characteristics  of  speckle  patterns  that  were 
observed  and  used  for  computer  model  verification  were  typi-  | 

cal  speckle  size  and  intensity  distribution  across  the  j 

speckle  pattern.  | 

i 

Average  speckle  size  is  directly  proportional  to  the  j 

laser  wavelength  and  range  and  inversely  proportional  to  the 
reflecting  aperture  size.  This  phenomenon  is  predicted  by 
theory  (see  theory  section) , was  shown  in  the  experimental 
investigation  of  speckle  patterns  (see  experimental  investi- 
gation section) , and  was  also  observed  in  the  calculations 
and  plots  from  the  speckle  model. 

The  intensity  distribution  in  the  far  field  diffraction 
pattern  is  very  complex  when  the  reflecting  surface  is  ran- 
domly rough.  To  verify  the  computer  model,  speckle  patterns  j 

were  first  calculated  for  a series  of  smooth  reflecting  sur- 
faces and  the  results  compared  very  well  with  theory.  For 
a square  complex  aperture  function,  the  first  minor  lobe 
amplitude  was  approximately  one  order  of  magnitude  lower 
than  the  central  maxima.  Each  succeeding  minor  lobe  was 
significantly  less  than  the  preceeding.  Vfhen  the  scattering 
surface  is  randomly  rough,  the  speckle  pattern  appears  as  a 
random  array  of  bright  and  dark  spots  and  minor  lobes  of  the 
diffraction  pattern  can  approach  the  amplitude  of  the  central 

I 

maxima.  This  phenomenon  was  observed  and  presented  in  the  I 
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experimental  section  of  this  thesis. 
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Thesis  Conclusion  Applicable  to  ALOT  System  Design 

The  last  chapter  of  this  thesis  gives  a summary  of  the 
thesis  project  results  that  are  applicable  to  the  ALOT  sys- 
tem design  study.  The  two  areas  that  are  summarized  are 
average  speckle  size  in  a far  field  diffraction  pattern  and 
the  spatial  frequencies  within  a speckle  pattern. 
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II. 


Diffraction  and  Fourier  Analysis  Theory 


Laser  light  differs  from  common  radiant  light  (sun, 
candle,  light  bulb,  etc.)  because  it  is  composed  of  a single 
frequency  and  displays  both  temporal  and  spatial  coherence. 
Temporal  coherence  expresses  the  length  of  time  that  the 
phase  of  the  wavefront  at  a given  point  in  space  can  be 
accurately  predicted.  Spatial  coherence  is  the  ability  to 
predict  the  phase  at  one  point  in  the  laser  beam  relative  to 
another  point,  either  longitudinally  or  laterally,  at  the 
same  instant  of  time.  For  a laser,  the  temporal  and  spatial 
coherence  characteristics  are  extremely  large  compared  to 
radiant  light. 

Finding  an  exact  explanation  and  mathematical  solution 
to  the  scattering  process  of  light  has  been  a major  area  of 
research  throughout  the  history  of  physics.  In  1678  the 
Dutch  physicist,  Christian  Huygens,  formulated  a theory  that 
stated:  (Ref. 7:60),  "every  point  on  a primary  wavefront 

serves  as  the  source  of  spherical  secondary  wavelets  such 
that  the  primary  wavefront  at  some  later  time  is  the  enve- 
lope of  these  wavelets.  Moreover,  the  wavelets  advance  with 
a speed  and  frequency  equal  to  that  of  the  primary  wave  at 
each  point  in  space."  In  the  early  1800  s Augustin  Jean 
Fresnel  recognized  that  the  huygens  principle  should  be 


3 


modified  to  include  the  wavelength  of  the  diffracted  light. 

He  corrected  this  apparent  deficiency  with  his  theory  which 
stated  (Ref. 7:330),  "every  unobstructed  point  of  a wavefront, 
at  a given  instant  of  time,  serves  as  a source  of  spherical 
secondary  wavelets  (of  the  same  frequency  as  the  primary 
wave) . The  amplitude  of  the  optical  field  at  any  point  be- 
yond is  the  superposition  of  all  these  wavelets  (considering 
their  amplitudes  and  relative  phases)."  This  statement  is 
now  referred  to  as  the  Huygens-Fresnel  principle  of  diffrac- 
tion. 

Although  the  phrases  of  the  Huygens-Fresnel  principle 
are  relatively  easy  to  understand,  the  mathematical  state- 
ment of  the  principle  is  very  complex.  Many  brilliant  theo- 
reticians and  experimentalists  such  as  Helmholtz,  Kirchhoff, 
Green,  Young,  Rayleigh,  Sommerfeld,  and  Fraunhofer  formu- 
lated reliable  approximations  enabling  the  complex  mathema- 
tical model  of  diffraction  to  be  reduced  to  a solvable 
equation  (Ref . 8 : 30-54) . 

Huygens-Fresnel  Equation 

Equation  1 is  one  of  many  mathematical  representations 
of  the  Huygens-Fresnel  diffraction  formula  (Ref. 8:60). 

+CO 

E<,(x,y)  = Exp(jkz)g  2)  }//e  (u,v) 

® ]Xz  2z  -00  a ' 


Exp{^(u^+v^)  }Exp{~j^^- (ux+vy)  }dudv 


(1) 


where : 


Electric  field  strength  in  observation 
plane  (v/m)  (complex  envelope) 

Electric  field  strength  in  scattering 
plane  (v/m)  (complex  envelope) 

Range  between  scattering  and  observation 
planes  (m) 

Laser  wavelength  (m) 

Propagation  number  = 2tt/X  (1/m) 

Equation  1 describes  the  relationship  between  the  complex 
amplitude,  Eg(x,y),  of  the  field  strength  of  the  scattered 
monochromatic  wave  in  the  observation  plane  due  to  a complex 
field  strength,  E^(u,v),  over  an  aperture  in  the  scattering 
plane.  See  Fig.  1 for  geometrical  layout.  Note  that  E^(u,v) 


Eg(x,y)  = 

Eg(u,v)  E 

z E 

A E 

k E 
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is  produced  by  the  monochromatic  wave  illuminating  the  scat- 
tering aperture.  It  will  be  shown  that  Eg.  1 will  simplify 
slightly  by  using  far  field  approximations  and  a complex 
aperture  function  model. 

Although  Eq.  1 was  formulated  through  the  use  of  exper- 
imental observations,  it  can  also  be  derived  from  the  differ- 
ential wave  equation.  An  equation  for  a scattered  field 
produced  by  a complex  aperture  function,  E^(u,v),  is  the 

a 

Fresnel-Kirchhof f diffraction  formula 

E (x,y)  = -i-//— i-JlxpljkrjlE  (u,v)dA  (2) 

® jA  A 

Equation  2 is  valid  for  small  0's  (see  Fig.  1)  in  the  range 
of  zero  to  15  degrees. 

An  approximation  of  r^  can  be  obtained  through  the  bi- 
nomial expansion  of 

r^  = (z^+ (x-u) (y-v) ^ ^ (3) 

When  this  expansion  is  completed,  only  those  terms  with 
powers  of  two  or  less  are  retained.  This  is  known  as  the 
Fresnel  approximation  and  is  valid  for 

(x-u)  ^+(y-v)  (4) 

^ MAX 
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with  the  previous  approximations,  the  Fresnel-Kirchhof f 
diffraction  formula  (Eq.  2)  reduces  to  the  Huygens-Fresnel 
diffraction  formula  (Eq.  1) . For  the  problems  investigated 
in  this  thesis,  Eq.  1 can  be  further  simplified  by  envoking 
the  Fraunhofer  far  field  approximation 


2 


MAX 


(5) 


This  approximation  insures  that  the  quadratic  phase  factor, 
from  Eq.  1, 


Exp{  jii(u^+v^  ) } (6) 

2z 

of  the  incident  radiation  is  approximately  unity  over  the 
aperture.  Goodman  states  (Ref. 8:61),  "Fraunhofer  diffraction 
patterns  can  be  observed  at  distances  closer  than  implied  by 
Eq.  5,  provided  the  aperture  is  illuminated  by  a spherical 
wave  converging  toward  the  observer,  or  if  a converging  lens 
is  properly  situated  between  the  observer  and  the  aperture." 

In  summary,  it  has  been  stated  (although  not  proven  in 
this  thesis)  that  the  Huygens-Fresnel  diffraction  formula 
(Eq.  1)  can  be  derived  from  the  differential  wave  equation 
if  realistic  approximations  are  applied  to  the  Fresnel- 
Kirchhoff  diffraction  formula  (Eq.  2) . In  addition,  Eq.  1 is 
further  simplified  by  the  Fraunhofer  far  field  approximation 
(Eq.  5)  to  provide  Eq.  7 on  the  following  page. 
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(7) 


Eg(x,y)  = E^xp  ( j_kz_)_  Exp{  lii(x^+y^  ) } 
j Xz  2z 

+ 00  _ ■ o 

/^E^(u,v)Exp{— i±-i(ux+vy)  jdudv 
X z 

Equation  7 represents  the  complex  amplitude  of  the  electric 
field  strength  (v/m)  distribution  at  an  observation  plane 
resulting  from  the  backscatter  of  an  electric  field  incident 
on  an  aperture.  All  the  quantities  in  Eq.  7 are  easily 
measured  or  defined  except  E^(u,v),  the  complex  aperture 
function,  when  the  aperture  is  a rough  reflecting  surface. 

Complex  Aperture  Function  Model 

The  complex  aperture  function,  Eg(u,v),  in  Eq.  7 is  a 
representation  of  the  complex  amplitude  of  the  electric  field, 
immediately  after  the  aperture,  that  propagates  through 
space  to  produce  a diffraction  pattern.  It  can  be  assumed 
that  Eg(u,v)  is  produced  by  an  infinite  number  of  fictitious 
spherical  point  sources  at  the  aperture,  where  each  source 
has  a unique  amplitude  and  phase.  Each  of  these  sources 
contributes,  by  addition  of  amplitudes  and  phases,  to  every 
point  on  the  far  field  diffraction  pattern.  Therefore,  the 
relative  phases  of  the  sources,  at  the  aperture,  directly 
effect  the  shape  of  the  diffraction  pattern.  The  funda- 
mental problem  then  is  to  construct  an  Ea(u,v)  function  that 
properly  describes  the  amplitude  and  phase  cf  every  point  in 
the  aperture. 

The  simplest  case  that  can  be  discussed  is  when  the 
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aperture  is  illuminated  by  a coherent  laser  beam.  The 
upper  drawing  in  Fig.  2 on  the  following  page,  illustrates 
the  geometry  of  the  problem.  The  complex  amplitude  of  each 
spherical  wavelet  can  be  described  by 

E^(u,v)  = (b)Exp{  j (kbj^)  } (8) 

From  the  geometry  shown  in  the  upper  drawing  of  Fig.  2 it  is 
obvious  that  at  b=0  all  relative  phases  are  zero  and  can  be 
assumed  equal  to  zero.  Thus,  Eq.  8 reduces  to 

E^(u,v)  = Eg  (b=0)  = Constant  (9) 

and  E^(u,v)  in  Eq.  7 can  now  be  moved  outside  the  double 
integral . 

When  the  aperture  becomes  a rough  reflecting  surface, 
as  in  the  operation  of  the  ALOT  pointing  system,  the  surface 
height  variations  can  cause  a randomness  in  the  phase  distri- 
butions of  Ea(u,v).  Rayleigh  (Ref. 9:9)  suggested  an  elemen- 
tary equation  to  determine  the  phase  at  each  point  of  the 
complex  aperture  function.  The  path  difference  (AC+CD-AB) 
between  waves  1 and  2 in  the  lower  drawing  of  Fig.  2 is 

A Path  = 2h(u,v)Cos3  (10) 

The  phase  difference  between  the  two  waves  is 

A0  = il(APath)  = 4^h(.u,_y)Cos^ 

X X 
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Equation  11  can  be  used  as  the  phase  modulation  of  the 
reflected  field.  To  complete  the  complex  aperture  function, 
E^(u,v),  terms  should  be  included  to  represent  the  amplitude 
of  the  field.  Let  P(u,v)  be  the  absolute  field  strength 
distribution  at  each  (u,v)  location  of  the  incident  radia- 
tion and  let  r be  the  reflectivity  of  the  scattering  surface. 
In  this  thesis  the  assumption  is  made  that  the  incident 
radiation  is  a plane  wave  from  a coherent  laser  where  the 
field  strength  is  a constant,  Eg  , across  the  beam  and  the 
reflectivity  is  a constant  across  the  aperture.  The  complex 
aperture  function  becomes 

E^(u,v)  = r Eg  Exp{i^h  (u,v)  CosB)  (12) 

Goodman  (Ref. 6: 65)  states,  "that  a simplistic  relation 
(for  Eg^(u,v))  often  used  in  analysis  and  reasonably  accurate 
if  the  surface  slopes  are  small."  is 

E^(u,v)  = r Eg  Exp{i^  (1+CosB)  h (u,v)  } (13) 

a A 

A quick  check  of  Eq.s  12  and  13  shows  that  the  phase  of  Eq.l3 
varies  more  slowly  than  that  of  Eq.  12,  however,  at  10 
degrees  or  less  inclination  (3)  of  the  incident  radiation 
they  differ  by  less  than  .033  radians.  Since  both  are  sim- 
plistic models  and  reasonably  equal  within  10  degrees  and 
Eq.  13  has  been  used  by  Goodman  in  the  study  of  laser 
speckle  statistics,  Eq.  13  will  be  used  in  this  thesis  as 
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the  model  of  the  complex  aperture  function,  E^(u,v). 

a 

with  the  addition  of  Eq.  13,  the  Fraunhofer  diffraction 
formula,  Eq.  7,  becomes 

E„(x,y)  = Exp{iii(x^+y^ ) }r  Eq 

^ jXz  2z 

4*00 

//Exp{l^(l+Cosg)  h (u,v)  }Exp{li^(ux+vy)  }dudv  (14) 

" X Az 

It  is  recognizable  that  the  double  integral  is  proportional 
to  the  two  dimensional  Fourier  transform  of  the  complex 
aperture  function,  Eg(u,v).  Before  any  additional  material 
is  covered  on  diffraction  patterns,  Fourier  analysis  theory 
of  discrete  data  processing  must  be  discussed. 

Fourier  Analysis 

This  subsection  of  theory  is  devoted  to  Fourier  analysis 
and,  in  general,  the  relationships  between  Fourier  inte- 
grals, discrete  Fourier  transforms  (DFT) , and  fast  Fourier 
transforms  (FFT) . In  addition,  the  provisions  for  the 
existence  of  each  transform  type  will  be  discussed  along 
with  the  pitfalls  associated  with  discrete  data  Fourier 
transforms.  The  DFT  is  defined  to  be 

N-1 

g(r)  =i  j^EqX  (k)  Exp{- j2iTrk/N}  (15) 
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where : 


g(r)  = coefficient  of  DFT  (sample) 

x(k)  = coefficient  of  spatial  series  (sample) 


N = Total  samples  in  spatial  series 
The  connection  between  discrete  and  continuous  Fourier  trans- 


form representations  are  shown  below. 


Discrete  Fourier  Transform 


The  Fourier  series  is  represented  by 


f(t)  = ^ + Z aj^Cos  (niTt/p) 
2 n=l 


+ Z bj^  Sin(mrt/p)  (16) 

n=l 


where 


^ 2p 

a,  = i / f(t)dt  (17) 

^ P 0 


1 

a = — / f (t)  Cos  (niTt/p)  dt  (18) 

u p 0 


1 

b_  = - / f (t)  Sin  (niTt/p)  dt  (19) 

“ p 0 


and  p is  the  period  of  the  time  series.  If  the  Dirichlet 
conditions  are  met  plus  the  integral  of  jf(t)  j from  minus  to 
plus  infinity  exists,  then  the  Fourier  series  (Eq.  16)  can 


i be  replaced  by  the  Fourier  integral.  Equations  20  and  21 

are  the  Fourier  integral  transform  pairs. 

I 

f 
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x(t)  = 7 g (f ) Exp{  j2iTf  t}df 


g(f)  = 7 X (t)  Exp{- j2irft}dt 


Equations  20  and  21  can  be  used  to  derive  the  DFT  (Ref. 10; 
280) . Start  with  Eq.  20  and  let 


t = t.  = iAt 
1 


i = 0,±1,±2, 


F = 1/At 


to  yield 


x(tj^)  = 7 g (f ) Exp{  j2TTif/F}df 


+CO  (k+l)F 

x{ti)  = , E f g (f ) Exp{  j2Trif/F}df 
k=-»  „ 


The  exponential  is  periodic  in  f with  period  F,  so  Eq.  25 
can  be  written  as 


x(iAt)  = f g (f)Exp{j27Tif/F}df 
0 P 


where 


g„(f)  = I g(f+kF) 
P k=-t» 


Since  gp(f)  is  periodic  with  period  F,  it  has  a complex 
Fourier  series  expansion  in  powers  of  the  exponential  where 
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1 


coefficients  are  given  by  (l/F)X(iAt),  as  shown  in  Eq,  26. 
Thus , 


g„(f)  = i ? X (iAt)  Exp{- j2iTif/F}  (28) 

f F i=— «= 


Now,  sample  in  the  frequency  domain  at  points  f =nAf, 
n=0,±l,±2, • • • , where  Af=l/T,  and  Eq.  28  becomes 


1 +0O 

g (nAf)  = ± Z x(iAt)Exp{-j2Trin/TF}  (29) 

^ F -i  =-“> 


If  TF=N  is  an  integer,  the  exponential  is  a periodic  func- 
tion of  i with  period  N.  Therefore 


1 N-1 

gp(nAf)  = i Xp  (iAt)  Exp{- j2irin/N}  (30) 


where 


+0O 

X (iAt)  = Z x((i+kN)At)  (31) 

P k=-«> 


The  sum  in  Eq.  31  gives  the  sampled  values  of 


X (t)  = , z”  x(t+KT)  (32) 

P k=-<» 


which  is  a periodic  function  of  t,  with  period  T,  and  is 

formed  from  x(t)  in  exactly  the  same  way  as  g (f)  is  formed 

P 

from  g(f)  (Eq.  21).  Substituting  At=T/N  for  1/F,  Eq.  30 

i becomes  Eq.  33  on  the  following  page.  Note  if  x(t)=0  for 

1 

I 

I 

\ 
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|tl>q,  Xp(t)=x(t)  (i.e.;  space  limited  function). 

N-1 

g (nAf)  = — E Tx^  (iAt)  Exp{- j2Trin/N}  (33) 

P N i=0  P 

which  is  the  DFT  (Eq.  15) . In  the  DFT  the  following  approx- 
imations apply: 


x(i)  = TXp(iAt) 

(34) 

g(n)  = gp(nAf) 

(35) 

This  proof  was  extracted  from  Rabiner  and  Rader 
and  provides  the  DFT  pairs  of 

(Ref .10:280) 

N-1 

g (n)  = - E X (k)Exp{-j2TTnk/N} 
P N k=0  P 

(36) 

N-1 

X (k)  = g (n)  Exp{  j2iTnk/N}  (37) 

P n=0  p 

and 

g (n)  = n^^  coefficient  of  DFT 
P 

X (k)  = k^^  coefficient  of  spatial  series 
P 

N = Total  samples  in  series 
k and  n = 0,  1,  2***N-1 

Equations  36  and  37  are  the  DFT  pairs  provided  Eqs.  20  and 
21  exist  and  the  approximations  of  Eqs.  34  and  35  are  used. 

The  FFT  is  an  efficient  method  for  calculating  the  DFT 
on  a computer.  This  thesis  will  not  discuss  the  development 
of  the  FFT  but  Ref.  10  part  2 gives  in  depth  details  of  the 
algorithm.  There  are  several  FFT  routines  that  are  used  in 
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the  signal  analysis  field,  however,  the  one  used  in  this 
thesis  and  most  widely  used  is  by  J.  W.  Cooley  and  J.  W. 
Tukey . 

The  number  of  multiplications  for  the  FFT  is  2NLog2N 
compared  to  for  the  DFT  for  a one  dimensional  data  array 
of  length  N samples.  Besides  this  computational  efficiency, 
the  FFT  uses  only  one  data  array  for  input  to  the  routine, 
internal  routine  calculations,  and  output  from  the  routine. 

There  are  two  major  problems  associated  with  using  the 
DFT  and  FFT  algorithms  compared  to  the  Fourier  integral; 
aliasing  and  leakage.  Aliasing  refers  to  the  fact  that 
high  frequency  components  of  a time  function  can  impersonate 
low  frequencies  if  the  sampling  rate  of  the  data  is  too  low 
and  therefore,  two  widely  different  functions  can  have  the 
same  DFT.  An  example  of  this  can  be  seen  in  Fig.  3 on  the 
following  page.  Curve  A is  the  time  function  sampled, 
however,  it  is  sampled  at  a frequency  lower  than  its  fre- 
quency. With  data  samples  available,  curve  B frequency  may 
be  assumed  since  two  or  more  data  samples  per  time  function 
period  are  required  to  determine  a frequency.  To  avoid 
aliasing  a sampling  rate  of  at  least  two  times  the  highest 
frequency  component  of  the  time  function  is  required.  This 
is  known  as  the  Nyquist  sampling  rate. 

When  working  with  systems  that  are  band-limited  it 
is  a good  practice  to  sample  at  more  than  twice  the  highest 
expected  frequency  of  the  function.  An  engineering  rule  of 
thumb  is  to  sample  at  a rate  of  four  to  ten  times  the 
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highest  expected  frequency.  The  maximum  sample  rate  is,  of 
course,  dependent  upon  the  data  handling  and  storage  charac- 
teristics of  the  overall  sampling  and  computer  system. 

The  second  problem  associated  with  DFT  is  that  of  leak- 
age resulting  from  a finite  length  record  of  data.  Consider 
the  continuous  sine  wave  shown  in  the  top  part  of  Fig.  4 on 
the  previous  page.  To  analyze  that  sine  wave,  it  is  sampled 
until  a finite  record  length  buffer  is  filled.  The  finite 
length  buffer  is  shown  in  the  middle  of  Fig.  4 and  the 
result  of  sampling  the  sine  wave  is  shown  in  the  bottom 
of  Fig.  4. 

In  the  frequency  domain,  the  Fourier  transform  of  the 
continuous  sine  wave  appears  as  shown  in  the  top  drawing  of 
Fig.  5 on  the  following  page.  This  diagram  shows  only  the 
positive  frequency  of  the  transform  and  it  is  an  impulse 
function  positioned  at  a frequency  equal  to  the  sine  wave 
frequency.  The  Fourier  transform  of  the  finite  length 
buffer  is  shown  in  the  middle  drawing  of  Fig.  5,  which  is 
a Sinc(x)  function.  The  result  of  sampling  the  sine  wave 
with  a finite  record  length  buffer  is  the  convolution,,  in 
the  frequency  domain,  of  the  impulse  function  and  the  Sine 
(x)  function,  as  shown  in  the  bottom  diagram  of  Fig.  5. 

This  result  is  known  as  windowing  and  produces  a maximum  peak 
at  the  frequency  of  the  sine  wave  and  also  a series  of  side- 
lobes  associated  with  the  finite  buffer  length. 

In  addition  to  the  problems  of  aliasing  and  leakage 
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associated  with  DFT  there  is  also  a frequency  resolution 
problem  with  the  data  returned  from  the  FFT.  The  top 
drawing  of  Fig.  6 on  the  following  page  represents  the  com- 
mon data  array  in  the  FFT  subroutine 

N = Maximum  array  dimension 
2q+l  E Total  non  zero  data  samples  in  array 
x(k)  = Magnitude  of  data  in  array  position  k 

where 


1 l-k-2q+l 

x(k)=  { 

0 2q+l<k-N  (38) 


r 


Recalling  the  DFT  equation  (Eq.  36)  and  replacing  the 
sequence  k = 0,1,2,***N-1  with  k = 1,2,3,***N,  the  result  is 


gp(n) 


L 

k=l 


(k)  Exp{  - j 27Tnk/N } 


(39) 


where  n=l,2,***N.  Since  x(k)  equals  zero  for  k>2q+l,  Eq.  39 
can  be  replaced  by 

2q+l 

g„(n)  = Z x^(k)  Exp{- j2irnk/N}  (40) 

^ k=l  F 

and  since  x(k)  = 1 for  all  l-k-2q+l,  as  shown  in  Fig.  6 on 
the  previous  page,  Eq.  40  is  approximated  by 


g(n)  = 2qSinc{  2TTnq/N}Exp{  - j 2Trn/N  ( l+q/2)  } (41) 


The  exponential  in  Eq.  41  results  from  the  pulse  not  cen- 
tered around  zero,  and  thus  is  a frequency  shift  in  the 
frequency  domain.  For  the  analysis  in  this  thesis,  the 
exponential  term  can  be  neglected. 

The  function  g(n)  of  Eq.  41  goes  to  zero  when  nq/N 
equals  an  integer.  Let  the  first  integer  be  one,  therefore 

nq/N  = 1 (42) 


or 


n = N/q 


(43) 
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Equation  43  shows  that  there  are  N/q  samples  between  0 and 
X,  as  shown  on  the  bottom  diagram  of  Fig.  6.  Therefore, 
to  provide  excellent  resolution  between  0 and  x on  the 
|sinc(x)  1^  function  returned  from  the  FFT,  N should  be  large 
and  q small.  As  an  example,  if  N = 128  and  q = 8,  there 
will  be  16  data  samples  between  0 and  the  first  null  on  the 
jsinc(x)  function. 

In  summary,  the  Fourier  analysis  review  has  shown  that 
the  DFT  and  FFT  are  excellent  representations  of  the  Fourier 
series  and  integrals  if  the  Dirichlet  conditions  are  met  and 
the  integral  of  ]f (t)  j exists.  To  accurately  represent  a 
continuous  time  function  with  discrete  sampling,  the  contin- 
uous function  should  be  sampled  at  a rate  greater  than  twice 
its  highest  frequency  component.  This  sampling  procedure 
will  reduce  or  eliminate  the  aliasing  problem.  In  addition 
finite  record  length  buffers  will  cause  leakage  and  corrupt 
the  data  returned  from  the  FFT.  Because  it  is  desirable  to 
have  a high  resolution  of  data  between  nulls  of  the  function 
returned  from  the  FFT,  the  actual  number  of  data  samples 
input  to  the  FFT  buffer  should  be  considerably  less  than  the 
total  size  of  the  input  buffer.  With  this  Fourier  transform 
review,  it  is  now  possible  to  complete  the  development  of 
the  Huygens-Fresnel  diffraction  formula  used  to  develop  the 
computer  model  of  speckle  in  this  thesis. 

Intensity  Distribution  at  an  Observation  Plane 

Equation  7 is  the  Huygens-Fresnel  diffraction  formula 


28 


I 

I 

i 


I 


1 


used  for  speckle  analysis  in  this  thesis.  This  form  of  the 
Huygens-Fresnel  equation  is  sometimes  called  the  Fraunhofer 
diffraction  equation  because  of  the  Fresnel  and  Fraunhofer 
assumptions  (Eqs.  4 and  5).  For  convenience,  Eq.  14  is 
repeated  in  Eq.  44. 

E (x,y)  = 5ii£li]i5l£xp{  jJl(x^+y^ ) }r  E 
jXz  2z  ° 

//Exp{J-^(l+Cos3)  h (u,v)  }Exp{~j-?-(ux+vy)  }dudv  (44) 

_oo  X Xz 


Although  not  evident,  Eq.  44  is  a time  varying  function 
and  the  intensity  distribution  can  be  obtained  from 
Eg (x ,y) Exp{ jwt}  by  taking  the  time  average.  Thus 

I (x,y)  = < (R^E(x,y) ) (45) 

where  < > denotes  time  average  which  reduces  to 

I(x,y)  = JjE  (x,y)  E*  (x,y)  (46) 

for  monochromatic  light. 

Since  Eq.  46  contains  the  electric  field  at  the  obser- 
vation plane  multiplied  by  its  complex  conjugate,  the  expo- 
nentials outside  the  integrals  of  Eq.  44  disappear.  The 
resulting  intensity  distribution  at  the  observation  plane  is 
given  by  Eq.  47  on  the  following  page,  where  I(x,y)  is  in 
WATTS/m^  and  is  sometimes  called  irradiance.  Recalling 
Fig.  1 shows  the  geometry  of  the  diffraction  problem. 
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I(x,y)  = h 1 — i{//Exp{2il(l+CosB)h(u,v) } 

X Z -00  X 

Exp{~j^^ (ux+vy) }}dudv| ^ (47) 

Xz 

The  proper  scaling  per  unit  at  the  observation  plane 
can  be  determined  by  solving  Eq.  47  for  a simplified  aper- 
ture function.  Let  the  aperture  be  square  with  linear 
dimension  b the  total  width.  Also  let  r and  h(u,v)  equal 
one  across  the  aperture  with  a B of  zero.  From  Eq.  47, 
I(x,y)  becomes 


I(x,y)/I^ 


fc-; ,TTbx.  I 

2 

TTbX 

TTby 

Xz 

Xz 

(48) 


The  distance  out  to  the  first  null  in  I(x,y)  is  determined 
by 


TTbx  ^ 
Xz 


IT 


or 


and  also 


(49) 


(50) 


(51) 
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Although  not  proven  here  (Ref. 7: 352)  the  intensity 
pattern  produced  by  a circular  aperture  is 


I(R)/1, 


(2iTrR/Xz)  ' 2 
irrR/Xz 


where 

r = Radius  of  aperture 
R H Radius  out  in  intensity  plane 
Jj  E Bessel  function  of  first  order 


(52) 


and  the  distance  out  to  the  first  null  of  the  intensity 
pattern  is  given  by 


R = ^ (53) 

2r 

From  Eq.  43  of  the  FFT  discussion,  there  are  n data 
values  between  zero  and  the  first  null  in  the  intensity  pat- 
tern returned  from  the  FFT  subroutine  (n  = N/q) . Recall 
N was  the  total  samples  in  the  input  buffer  and  q was  the 
total  samples  of  real  data  in  that  buffer.  Therefore,  the 
linear  distance  per  sample  at  the  observation  plane  is 

Ax  = Ay  = x/n  = y/n  = (54) 

N/q 

However,  b was  the  total  pulse  width  subdivided  into  q 
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samples,  therefore 


Ax  =Ay  = Al  . a (55) 

N b 

Equation  55  gives  the  linear  distance  per  sample  of  the 
speckle  pattern  in  the  observation  plane  when  the  input  is 
a square  aperture.  It  is  obvious  that  the  resolution  dis- 
tance at  the  observation  plane  is  directly  proportional  to 
A and  z and  inversely  proportional  to  N (total  record 
length  of  data  buffer) . Equation  47  will  be  programmed  to 
calculate  the  speckle  phenomenon  for  rough  surface  scatter- 
ing and  Eq.  55  to  determine  the  unit  sample  distance  at  the 
observation  plane. 

Summary  of  Theory  Section 

As  revealed  in  the  introductory  paragraphs  of  this 
chapter,  the  search  for  a mathematical  solution  to  the  dif- 
fraction of  light  began  300  years  ago.  In  the  early  1800's 
the  Huygens-Fresnel  equation  (Eq.  1)  was  developed  as 
hypothesized  from  experimental  data  and  later  proven  to  be 
a valid  solution  of  the  scalar  wave  equation.  The  general 
Huygens-Fresnel  equation  was  simplified  by  using  the  Fresnel 
approximation  to  replace  the  spherical  waves  by  plane  waves. 
Equation  4 gives  the  conditions  when  the  approximation  is 
valid.  The  Fraunhofer  far  field  approximation  (Eq.  5)  was 
also  utilized  to  remove  the  quadratic  phase  factors  at  the 
aperture.  The  generalized  Huygens-Fresnel  equation  (also 
called  Fraunhofer  diffraction  equation)  was  then  revised  to 
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include  rough  surface  scattering.  The  term  added  to  the 
equation  is  called  the  complex  aperture  function  (Eq.  13) 
as  suggested  by  Goodman  (Ref. 6).  The  final  result  was  an 
equation  that  predicted  the  electric  field  strength  (v/m)  at 
an  observation  plane  when  an  incident  laser  beam  was  scat- 
tered from  a rough  surface  (Eq.  14)  . 

Before  the  intensity  field  (irradiance)  formula  was 
developed,  a detailed  investigation  of  Fourier  analysis  was 
accomplished.  It  was  shown  that  the  Fourier  integrals, 
discrete  Fourier  transforms,  and  fast  Fourier  transforms 
are  analogous  if  several  conditions  exist.  As  stated,  to 
represent  a time  function  by  the  Fourier  integral  instead 
of  the  series,  the  Dirichlet  conditions  must  be  met  and  the 
integral  of  |f (t) | must  exist  for  all  time.  The  discrete 
and  fast  Fourier  transforms  are  valid  approximations  to  the 
integral  transform,  however,  the  time  function  must  be 
sampled  at  the  Nyquist  rate  or  greater  to  remove  the  alias- 
ing problem.  Lea)cage  was  covered  in  detail  and  it  was  shown 
that  discrete  and  fast  Fourier  transformed  data  are  cor- 
rupted by  the  use  of  finite  record  length  data  buffers.  The 
Fourier  analysis  section  was  concluded  by  proving  that  the 
input  array  to  the  fast  Fourier  transform  should  be  large 
and  the  actual  data  samples  should  fill  only  a portion  of  the 
array.  This  restriction  was  required  to  provide  high  reso- 
lution of  the  I Sine (x)  envelope  between  zero  and  the  first 
null  at  the  observation  plane. 

The  final  section  of  the  theory  chapter  developed  the 
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equation  to  predict  the  intensity  pattern  (Eq.  47)  at  an 

i 

observation  plane.  As  stated,  the  intensity  pattern  is  | 

actually  proportional  to  the  magnitude  squared  of  the  two  ; 

j 

dimensional  Fourier  transform  of  the  complex  aperture  func-  ; 

tion.  The  last  equation  to  be  developed  (Eq.  55)  provided  i 

the  proper  scaling  of  distance  at  the  observation  plane. 

Equation  47,  the  far  field  diffraction  equation,  was  used  j 

to  develop  the  computer  model  for  the  speckle  phenomenon.  j 


i 
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III.  Experimental  Investigation  of  Speckle  Patterns 


T 


A laboratory  experiment  was  conducted  to  measure  the 
intensity  distribution  in  speckle  patterns  and  identify  the 
trends,  from  pattern  to  pattern,  as  the  incident  laser  beam 
size  and  range  between  scattering  and  observation  planes 
were  changed.  By  observing  and  cabulating  these  quantities 
and  then  comparing  them  to  the  theoretical  computer  model, 
a qualitative  estimate  of  the  model  "goodness"  was  estab- 
lished. This  section  of  the  thesis  covers  only  the  experi- 
mental investigation,  while  the  comparison  of  the  computer 
model  patterns  with  the  experimental  data  is  covered  in  a 
later  chapter. 

Experimental  Set-Up 

Figure  7 shows  the  layout  of  the  laboratory  equipment 
used  in  the  experimental  investigation  of  the  speckle  pat- 
terns. The  equipment  used  was  available  in  the  physics 
laboratories  of  the  Air  Force  Institute  of  Technology  at 
Wright-Patterson  Air  Force  Base,  Ohio.  Table  I,  lists  the 
equipment  shown  in  Fig.  7.  It  should  be  noted  that  measure- 
ments of  the  reflecting  surface  roughness  were  not  obtained 
and  the  consequences  of  this  are  discussed  later  in  this 
chapter. 

Experimental  Procedure 

To  provide  an  understanding  of  the  diffraction  proper- 
ties of  speckle  patterns  two  experimental  procedures  were 
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Fig.  7 Laboratory  Equipment  Layout 


Table  I 


Laboratory  Equipment  Nomenclature 

Laser:  HeNe  50  mW  spectra  physics  model  125  ! 

X = 6328  8 ; 

Laser  collimator:  Model  336  multiple  wave- 

length collimating  lens,  spectra 
physics 

Reflecting  surface:  Oxidized  aluminum  scrap 

metal  (7.5  x 7.5  x 0.5  cmj 

Photomultiplier:  Photomultiplier  micropho- 
tometer No. 10-213 

Photomultiplier  tube:  Type  "S",  EMI 

Focusing  lens:  4 cm  diameter,  20  cm  focal 

length 

Traversing  micrometers:  Standard  laboratory 

traversing  micrometer  with  1 ym. 
resolution 

Camera:  Graphix,  aperture  and  focal  plane 

shutters 

Film:  Polaroid,  black  and  white,  Polapan/ 

type  52 

developed.  The  first  procedure  was  to  observe  the  speckle 
patterns  at  various  observation  planes  while  keeping  the 
beam  focused  on  the  reflecting  surface  (smallest  incident 
beam  spot) . Features  of  these  speckle  patterns  (a  pattern 
formed  by  a series  of  bright  and  dark  spots)  were  identi- 
fied and  could  be  recognized  at  several  observation  planes. 

Photographs  of  one  section  of  each  speckle  pattern,  at 
three  observation  planes,  were  obtained.  For  the  next  part 
of  the  experiment  the  scattering  plane  was  traversed  paral- 
lel to  the  incident  radiation.  This  operation  defocused 
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the  radiation  on  the  scattering  surface  which  increased  the 
spot  size,  however,  a constant  distance  between  the  scat- 
tering and  observation  planes  was  maintained.  Three  scat- 
tering plane  locations  were  used  and  a photograph  of  a 
section  of  each  speckle  pattern  was  obtained. 

The  second  procedure  was  similar  to  the  first  except  a 
photomultiplier  was  used  to  measure  the  speckle  pattern 
intensities  (irradiance) . Use  of  the  photomultiplier  pre- 
sented some  problems  that  had  to  be  solved.  The  first 
problem  was  the  enormous  area  of  the  input  aperture  (approx- 
imately 25  sq.  cm.)  of  the  photomultiplier.  Since  this 
area  encompassed  many  bright  and  dark  spots  of  the  speckle 
pattern,  speckle  intensity  (actually  a measurement  of  the 
electric  field  strength)  resolution  was  very  poor  and,  in 
fact,  an  average  was  obtained.  This  problem  was  solved  by 
decreasing  the  aperture  of  the  photomultiplier  tube  with  a 
piece  of  metal  that  contained  a 2.50  ± 0.05  mm.  diameter 
hole.  Tin  foil  was  then  taped  over  the  metal  and  a pin 
hole  (0.75  ± 0.05  mm.  diameter)  was  made  to  coincide  with 
the  hole  in  the  metal  plate.  Because  of  the  reduction  in 
the  aperture  area  of  the  photomultiplier,  intensity  measure- 
ments obtained  were  relative  and  not  calibrated.  The  second 
major  problem  with  using  the  photomultiplier  to  measure 
relative  intensity  strengths  of  the  speckle  patterns  was  the 
large  amount  of  data  samples  required  to  characterize  a 
pattern.  The  solution  to  this  problem  was  to  search  the 
speckle  pattern,  at  the  observation  plane,  for  the  maximum 
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Speckle  patterns  were  compared  when  the  range  was  held 
constant  and  the  reflecting  aperture  size  was  changed.  The 
smallest  incident  radiation  spot  size  was  obtained  when  the 
scattering  surface  was  placed  at  the  focal  plane  of  the 
focusing  lens.  The  spot  size  increased  when  the  scattering 
surface  was  moved  5 mm.  from  focus  condition  and  also  at 


10  mm.  off  focus.  In  all  cases,  the  incident  radiation 
spot  size  was  not  accurately  measured.  The  object  of  this 
part  of  the  experiment  was  to  observe  the  trend  in  the 
speckle  size  as  the  aperture  size  increased. 


Photographic  Evaluation  of  Speckle  Patterns 

Photographs  are  an  easy  way  to  record  speckle  patterns 
in  the  laboratory,  however,  measurement  of  contrast  is 
usually  very  hard  to  accomplish  from  photographs.  Figure  8 
shows  a speckle  pattern  photographed  0.5  ± 0.005  m.  from  the 
scattering  plane  with  the  incident  laser  beam  radiation 
focused  on  the  scattering  surface.  Areas  of  intensity 


39 


Fig.  8 

Speckle,  rough 


surface  scattering, 
range=0,5  m 


Fig.  9 

Speckle,  rough 
surface  scattering, 
range=1.0  m 


Fig.  10 

Speckle,  rough 
surface  scattering, 
range=l . 5 m 
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Fig . 

Speckle,  rough 
surface  scattering, 
range=0.5  m,  aper- 
ture less  than 
0.^  mm  diameter 


Fig.  12 

Speckle,  rough 
surface  scattering, 
range=0.5  m,  aper- 
ture less  than 
0.5  mm  diameter 


Fig.  13 

Speckle,  rough 
surface  scattering, 
range=0.5  m,  aper- 
ture less  than 
0 . 8 mm  diameter 
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enhancement  and  cancellation  are  readily  apparent,  and  in 
addition,  areas  with  unique  speckle  shapes  can  be  located. 

The  centimeter  scale  overlay  on  Fig.  8 highlights  a speckle 
and  shows  its  length  to  be  14.5  ± 0.5  mm.  Figure  9 is  a 
photograph  of  the  same  speckle  pattern  area  with  the  obser- 
vation plane  at  1.0  + 0.005  m.  The  speckle  is  now  29.5 
± 0.5  mm.  long.  Figure  10  is  a photograph  of  the  same  area 
at  1.5  + 0.005  m.  with  the  speckle  now  43.5  ± 0.5  mm.  long. 

These  three  photographs  illustrate  one  of  the  prin- 
ciples of  diffraction  theory.  Recalling  from  the  theory 
section,  (Eqs.  52  and  53)  the  linear  extent  of  a typical 
speckle  lobe  is  directly  proportional  to  the  range  between 
the  scattering  and  observation  planes.  The  measurement 
values  from  the  three  figures  support  the  diffraction  theory, 
i.e.,  when  the  range  is  doubled  or  tripled  so  does  the 
speckle  spot  diameter. 

For  the  next  experimental  investigation,  the  range  from 
observation  plane  to  scattering  plane  was  held  constant  at 
0.5  ± 0.005  m.  Figure  11  is  a photograph  of  a section  of 
the  speckle  pattern  when  the  incident  laser  beam  radiation 
was  focused  on  the  scattering  surface.  Although  an  accurate 
measurement  of  the  scattering  spot  size  was  not  obtained, 
it  was  less  than  0.2  mm.  diameter.  Figures  12  and  13  are 
photographs  of  the  same  speckle  pattern  area  as  in  Fig.  11 
with  the  scattering  plane  moved  5.0  and  10.0  ± 0.5  mm.  from 
the  focus  position.  The  incident  radiation  diameters  at 
these  locations  were  less  than  0.5  and  0.8  mm.  respectively. 
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It  is  very  difficult  to  make  any  quantitative  statement 
about  the  characteristics  shown  in  Figs.  11  through  13, 
however,  qualitatively  the  speckle  sizes  were  reduced  as  the 
aperture  size  increased.  To  obtain  more  quantitative  infor- 
mation about  the  speckle  phenomenon,  increased  resolution  of 
intensity  contrast  in  the  speckle  pattern  was  required.  A 
photomultiplier  was  used  to  provide  the  increase  in  reso- 
lution. 

Photomultiplier  Evaluation  of  Speckle  Patterns 

As  shown  in  Fig.  7,  the  photomultiplier  entrance  (pin 
hole)  aperture  was  positioned  at  the  speckle  pattern  obser- 
vation plane.  After  a maximum  intensity  in  the  speckle 
pattern  was  located,  the  photomultiplier  apparatus  was 
aligned  to  traverse  parallel  to  the  scattering  surface 
through  the  maximum  intensity.  As  stated  earlier,  the  photo- 
multiplier aperture  position  was  incremented  in  0.5  mm. 
steps.  After  each  increment,  an  indicated  voltage  (propor- 
tional to  the  intensity  of  the  speckle  pattern  over  by  pin 
hole  was  recorded  in  the  laboratory  notebook.  This  proce- 
dure was  repeated  for  four  speckle  patterns;  1.)  Range  one 
meter,  scattering  surface  at  focus  plane,  2.)  Range  two 
meters,  scattering  surface  at  focus  plane,  3.)  Range  two 
meters,  scattering  surface  five  millimeters  from  focus  plane, 
and  4.)  Range  two  meters,  scattering  surface  ten  millimeters 
from  focus  plane. 

To  aid  in  surveying  the  intensity  data,  a computer 
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program  was  written  to  normalize  the  data  and  provide  plots. 
Figures  14  through  17  show  the  normalized  speckle  patterns 
for  the  four  cases  described  in  the  previous  paragraph.  To 
obtain  the  frequency  content  of  each  speckle  pattern  a pro- 
gram was  written  to  perform  the  Fourier  transform  on  each 
recorded  speckle  pattern  and  to  produce  a frequency  plot  for 
analysis.  Figures  18  through  21  show  the  Fourier  trans- 
formed data  for  the  four  speckle  patterns  recorded. 

The  speckle  patterns  were  sampled  every  0.5  mm.  or  at 
a rate  of  2 cycles/mm.  Recalling  from  the  theory  section 
that  a sample  rate  of  twice  the  highest  frequency  content 
of  the  data  is  required  to  prevent  aliasing,  the  sample  rate 
used  provided  valid  frequency  information  below  1 cycle/mm. 
Inspection  of  the  Fourier  transform  plots  (Figs.  18  through 
21)  reveal  that  the  frequency  content  of  the  speckle  pat- 
terns was  far  below  the  Nyquist  rate.  Thus,  the  sample 
rate  was  much  greater  than  required  to  characterize  the 
speckle  patterns. 

Analysis  of  the  Speckle  Patterns 

A FFT  routine  was  used  to  calculate  the  Fourier  trans- 
form of  the  speckle  pattern  data  collected  with  the  photo- 
multiplier. The  Fourier  transform  of  this  intensity  data 
provided  the  spatial  frequency  content  of  the  segment  of 
each  speckle  pattern  sampled.  Although  each  speckle  pattern 
segment  was  only  a small  portion  of  the  total  speckle 
pattern,  the  spatial  frequencies  obtained  should  be 
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EXPERIMENTAL  INTENSITY  DATA 


Photomultiplier  Data  of  Speckle  Pattern 


Photomultiplier  Data  of  Speckle  Pattern 


EXPERIIiENTRL  INTENSITY  DflTR 


Fig.  16  Photomultiplier  Data  of  Speckle  Pattern 
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PRRnUFlNCY  CONTENT  OP  SPF.CKLE  PATTERN  RT  13.9  CYaES/hETEf^.UNIT 

Fig.  19  Fourier  Transform  of  Fig.  15  Intensity  Data 
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representative  of  the  total  pattern. 

To  explain  the  analysis  procedure.  Fig.  17  (intensity 
data)  and  Fig.  21  (Fourier  transform  of  the  data  on  Fig.  17) 
will  be  used.  Note  that  on  the  Fourier  transform  plot, 

Fig.  21,  there  is  a high  spatial  frequency  content  around 
zero  cycles/m.  This  amplitude  at  zero  frequency  results 
from  the  Fourier  transform  of  the  zero  frequency  level  of  the 
input  intensity  data.  Fig.  17.  The  three  REF  marks  on  the 
spatial  frequency  plot.  Fig.  21,  highlight  three  dominant 
frequencies  that  were  present  in  the  segment  of  the  speckle 
pattern  sampled.  Note  that  there  was  a frequency  peak  at 
± 37  units  on  Fig.  21  that  was  not  highlighted.  That  fre- 
quency was  not  highlighted  because  its  normalized  amplitude 
was  less  than  0.01,  and  only  peaks  greater  than  0.01  were 
referenced.  Since  the  amplitude  and  spread  of  a frequency 
peak,  on  the  Fourier  transform  plot,  are  a measure  of  the 
dominance  of  that  frequency  in  the  spatial  data,  the  0.01 
value  is  an  analyst  benchmark.  Since  the  Fourier  transform 
plots  do  not  provide  the  resolution  necessary  to  pick  the 
frequency  amplitudes  which  are  greater  than  0.01,  the  tabu- 
lated frequency/amplitude  data  from  the  program  was  used. 

A typical  speckle  width  was  also  calculated  for  each 
speckle  pattern  segment  analyzed.  To  obtain  the  typical 
speckle  width,  Eq.  56  was  used. 

Typical  speckle  width  = i (56) 

Highest  REF  frequency 
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since  the  highest  REF  frequency  was  an  analyst  choice  based 
on  normalized  amplitudes,  the  typical  speckle  width  may  be 
subject  to  error  and  should  be  verified.  As  an  example  use 
Figs.  16  and  20.  The  REF  frequency  on  Fig.  20  is  16  units 
or  222.4  cycles/m.  Using  Eq.  56  this  gives  a typical 
speckle  width  of  4.5  ±0.1  mm.  or  9.0  ± 0.2  units  on  Fig.  16. 
A close  inspection  of  the  intensity  data  shown  on  Fig.  16 
reveals  that  9 units  is  a good  typical  speckle  width  for  the 
data.  Figures  17  and  21  also  reveal  that  Eq.  56  is  a good 
evaluation  of  the  intensity  data.  Table  II  summarizes  the 
photomultiplier  data  analysis. 


Discussion  of  Experimenual  Results 

The  most  tangible  conclusion  from  the  experimental 
investigation  of  speckle  patterns  was  that  when  the  range 
from  scattering  to  observation  planes  was  doubled  or  tripled 
so  did  the  typical  speckle  size.  This  phenomenon  was  sup- 
ported by  theory  (Eqs.  52  and  53),  Figs.  8,  9,  10,  and  the 
results  shown  in  Table  II.  Theory  also  predicts  that  speckle 
size  should  decrease  as  the  reflecting  aperture  size  in- 
creases. Although  accurate  measurements  of  the  scattering 
spot  sizes  were  not  obtained,  the  experimental  results  sup- 
port the  theory.  This  phenomenon  can  be  seen  in  Figs.  11, 

12,  13,  and  the  results  shown  in  Table  II. 

Although  the  surface  roughness  of  the  scattering  pZane 
was  not  measured,  it  was  obvious  that  the  height  variations 
were  greater  than  the  laser  wavelength.  Recalling  from  the 
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theory  section,  when  the  reflecting  aperture  is  smooth,  the 
far  field  diffraction  pattern  is  an  Airy  pattern  where  the 
minor  lobes  are  more  than  an  order  of  magnitude  below  the 
central  major  lobe.  When  the  scattering  surface  is  rough, 
height  variations  greater  than  the  laser  wavelength,  the 
secondary  lobes  can  approach  the  amplitude  of  the  central 
peak.  This  characteristic  is  observed  in  Fig.  16  where  it 
is  very  easy  to  pick  out  the  diffraction  lobes.  From  Eq.  53, 
of  the  theory  section,  the  linear  dimension  of  a speckle 
produced  by  a circular  aperture  is 


R = 1- 22Xz 
2r 


(57) 


or  the  diameter  of  the  aperture  is 

D = 2r  = (58) 

R 

From  Table  II,  experiment  3,  if  we  use  the  typical  speckle 
size  (R)  of  4.5  mm.,  range  (z)  of  2m.,  and  wavelength  (X) 
of  6328  8,  and  Eq.  58,  the  aperture  diameter  should  Le 
approximately  0.3  mm.  Note  that  this  value  is  in  the  range 
of  the  estimated  spot  size.  Rough  surface  diffraction  will 
be  discussed  in  greater  detail  in  the  computer  model  veri- 
fication chapter. 


IV,  SPECKLE  Program 


A principle  part  of  the  time  allotted  to  the  completion 
of  this  thesis  was  devoted  toward  the  development  of  a reli- 
able and  accurate  computer  code  to  calculate  and  display 
speckle  (diffraction)  patterns.  This  program  was  written 
to  satify  two  requirements:  1,)  to  provide  a simulation 

tool  for  the  systems  control  engineers  involved  in  the  ALOT 
system  investigation  and  2.)  to  provide  visual  aids  on 
Fourier  optics  and  diffraction  theory  for  the  electro-optics 
systems  engineer. 

The  SPECKLE  program  was  designed  to  calculate  and  plot 
the  far  field  diffraction  pattern  resulting  from  the  rough 
surface  scattering  of  a laser  beam.  The  user  of  this  pro- 
gram must  supply  two  subroutines.  The  first  subroutine, 
HEIGHS,  provides  the  SPECKLE  program  with  input  data  includ- 
ing an  array,  DATA,  that  contains  the  scattering  surface 
height  measurements  normalized  by  the  laser  wavelength.  The 
arrangement  of  the  normalized  height  data,  in  the  array 
DATA,  defines  the  aperture  cross  section.  The  second  sub- 
routine, COMAPE,  calculates  the  complex  aperture  function 
and  returns  the  result  to  the  SPECKLE  program  through  the 
array  DATA.  In  the  COMAPE  subroutine,  the  normalized  sur- 
face height  data  is  modified  by  the  complex  aperture  func- 
tion model  which  includes  the  incident  radiation  field 
distribution,  aIFIEL  (u,v) , and  the  surface  reflection  coef- 
ficient, REF  (u,v) . These  subroutines  along  with  the 


57 


SPECKLE  program  are  explained  in  detail  in  this  chapter 

Overview  of  SPECKLE  Program 

The  SPECKLE  program  was  coded  in  Extended  Fortran  and 
designed  to  run  on  the  CD  6600  computer  at  Wright-Patterson 
AFB.  Dayton,  Ohio.  The  core  requirements  are  large,  245,000 
(octal)  words,  with  the  following  approximate  assignment  for 
major  routines: 


SPECKLE 

DATA 

FOURT 

DISSPLA 


Main  program 
Array 

FFT  subroutine 
Plotting  package 


75,000g  words 
100,000s  words 
2,000g  words 
67,000s  words 


The  DATA  array  is  used  for  input,  common  storage  for  FOURT 
and  DISSPLA,  and  output  routines.  This  is  a two  dimensional 
complex  array,  thus  requiring  two  consecutive  core  loca- 
tions (real  and  imaginary  parts)  for  every  data  value.  The 
overall  dimensions  of  this  array  determines  the  frequency 
resolution  of  the  data  returned  from  the  FOURT  subroutine 
(see  theory  section  of  this  thesis) . The  DISSPLA  package 
is  a library  program  available  on  the  CD  6600  computer  at 
Wright-Patterson  AFB. 

The  computer  code  consists  of  a main  program  (SPECKLE) 
and  six  subroutines  which  are  listed  on  the  following 
page . 
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r 1 


HEIGHS 

User  supplied  subroutine 
surface  height  data 

to 

input 

COMAP E 

User  supplied  subroutine  to 
culate  the  complex  aperture 

cal- 

function 

FOURT 

FFT  subroutine 

SHIFT 

Moves  data  in  array 

DATA 

DISSPLA 

Plotting  package 

VALU 

Function  subroutine 
DISSPLA 

to  assist 

The  following  paragraphs  describe  the  operations  of  the  sub- 
routines, the  input  data  descriptions,  and  finally  the  flow 
chart  and  operation  of  the  main  program. 

Subroutine  HEIGHS 

Subroutine  HEIGHS  is  a user  supplied  routine  that  in- 
puts data  to  the  SPECKLE  program.  The  following  is  a list 
of  variables  that  must  be  supplied  by  HEIGHS: 


; DATA  (IROW 

, ICOL)  Normalized  surface  height 
measurements 

IROW 

■ 

Maximum  row  size  of  DATA  array 

. ICOL 

Maximum  column  size  of  DATA  array 

MRDAT 

Maximum  number  of  height  data 
samples  in  u direction  (row) 

MCDAT 

Maximum  number  of  height  data 
samples  in  v direction  (column) 

ALAMBA 

Laser  wavelength  in  meters 

i RANGE 

Distance  from  scattering  to  obser- 
vation plane  in  meters 

i 
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AMESC  Linear  distance  between  surface 

height  samples  in  meters 

IN2D  Set  to  1 if  a two  dimensional  plot 

of  surface  height  data  is  desired, 
0 otherwise 

IN3D  Set  to  1 if  a three  dimensional 

plot  of  surface  height  data  is 
desired,  0 otherwise 

INT2D  Set  to  1 if  a two  dimensional  plot 

of  speckle  pattern  is  desired, 

0 otherwise 


INT3D  Set  to  1 if  a three  dimensional 

plot  of  speckle  pattern  is 
desired,  0 otherwise 


ISFREQ  Set  to  1 if  a plot  of  the  Fourier 

transform  of  the  two  dimensional 
speckle  pattern  is  desired, 

0 otherwise 


ISQAP  1 signifies  a 

ICRAP  1 signifies  a 

shape 

ICOAP  1 signifies  a 

shape 

ISMSU  1 signifies  a 

surface 

IROSU  1 signifies  a 

surface 


square  aperture  shape 
circular  aperture 

complex  aperture 

smooth  scattering 

rough  scattering 


Subroutine  HEIGHS  is  called  from  the  SPECKLE  program  by 


CALL  HEIGHS  (IROW,  ICOL,  MRDAT,  MCDAT,  ALAMBA, 
RANGE,  AMESC,  IN2D,  IN3D,  INT2D,  INT3D, 
ISFREQ,  ISQAP,  ICRAP,  ICOAP,  ISMSU,  IROSU) 


The  array  DATA  is  passed  through  common  by 
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COMPLEX  DATA 

COMMON  /A/  DATA  (128,  128) 

The  information  that  is  transferred  to  the  SPECKLE 
program  through  DATA  is  the  normalized  surface  height  mea- 
surements, i.e.,  actual  surface  height  measurements  divided 
by  the  laser  wavelength.  In  addition,  the  normalized  height 
data  is  centered  in  the  DATA  array  and  the  data  arrangement 
describes  the  aperture  cross  section.  The  array  DATA  is 
complex,  however,  the  surface  height  measurements  are  real, 
therefore,  the  data  must  be  entered  into  the  real  part  of 
each  DATA  location.  As  an  example,  if  the  size  of  the  DATA 
array  was  128x128  and  a square  aperture  cross  section  15x15 
samples  was  to  be  represented,  only  the  central  15x15  sam- 
ples of  the  DATA  array  would  contain  normalized  surface 
height  data.  All  other  locations  in  DATA  will  be  filled 
with  zeros. 

The  maximum  size  of  DATA  must  be  specified  by  IROW  and 
ICOL  at  128x128,  64x64,  32x32,  or  16x16.  The  maximum  number 
of  data  samples  across  the  aperture  cross  section  is  provid- 
ed to  the  main  program  by  MRDAT  and  MCDAT.  MRDAT  and  MCDAT 
must  be  less  than  IROW  and  ICOL,  and  to  provide  good  reso- 
lution in  the  data  returned  from  the  fast  Fourier  transform, 
the  ratios  of  IROW  to  MRDAT  and  ICOL  to  MCDAT  should  be 
large.  As  an  example,  let  IROW  and  ICOL  equal  128  and  MRDAT 
and  MCDAT  equal  15  or  31. 
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Subroutine  COMAPE 


COMAPE  is  a user  supplied  subroutine  that  transforms  the 
normalized  surface  height  data  into  the  complex  aperture 
function  distribution 

ADATA(I,J)  = REP(I,J) *AIFIEL(I,J)  * 

Exp{  j2iT*DATA{I,  J)  * (l.-Cos  (GLAZ)  ) } 

ADATA(I,J)  =0  if  DATA(I,J)  is  zero  (59) 

The  complex  aperture  function  distribution  is  returned  to 
the  SPECKLE  program  through  the  array  DATA(I,J).  The  input 
variables  to  the  complex  aperture  function  are  as  follows: 

REF(I,J)  Surface  reflection  coefficient 

distribution 

AIFIEL(I,J)  Incident  field  radiation  distri- 

bution (v/m) 

GLAZ  Angle  between  incident  field 

normal  and  scattering  surface 
normal  (-0.26  radians) 

The  variables  REF  and  AIFIEL  permits  the  user  to  vary 
the  reflection  coefficient  across  the  aperture  and  provide 
an  incident  field  profile  other  than  a plane  wave.  For  the 
examples  seen  in  this  thesis,  REF(I,J)  and  AIFIEL(I,J)  were 
set  to  constant,  REF  and  AIFIEL,  to  reduce  the  core  require- 
ments . 

Subroutine  SHIFT 

Subroutine  SHIFT  was  designed  to  rearrange  data  in  the 
complex  two  dimensional  array  DATA.  The  data  returned  from 
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the  FOURT  subroutine  (FFT)  has  the  zero  frequency  magnitude 
of  the  Fourier  transform  positioned  in  DATA  (1,1).  Recall- 
ing that  IROW  and  ICOL  define  the  maximum  dimensions  of  the 
array  DATA,  the  highest  frequency  of  the  transform  is  posi- 
tioned in  DATA  (1,  ICOL/2)  in  one  direction  and  DATA 
(IROW/2,  1)  in  the  other  direction.  The  remaining  three- 
fourths  of  the  array  DATA  contain  redundant  information. 

SHIFT  rearranges  the  data  to  position  the  zero  frequency 
component  in  the  center  of  the  array  and  provide  a symmet- 
rical Fourier  transform  function  around  DATA  (IROW/2 , ICOL/2 ) . 

Subroutine  FOURT 

The  FOURT  subroutine  (multi-dimensional  FFT)  used  in 
this  program  was  obtained  from  the  Electrical  Engineering 
Department  of  the  Air  Force  Institute  of  Technology,  Wright- 
Patterson  AFB,  Ohio.  The  subroutine  performs  the  Cooley- 
Tulcey  fast  Fourier  transform 

IROW  ICOL 

TRANSFORM{DATA(x,y) } = E E DATA(u,v) 

u=l  v=l 

Exp{-j27r{iiizili2lil  + (Y-l).(v-l,).}}  (60) 

IROW  ICOL 

where  IROW  and  ICOL  are  the  maximum  row  and  column  dimen- 
sions of  the  array  DATA. 

Subroutine  DISSPLA 

The  DISSPLA  subroutine  package  is  present  on  the  CD6600 
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computer  user  library  at  Wright-Patterson  AFB,  Ohio.  This 
display  package  was  used  because  it  contained  the  three 
dimensional  plotting  routines  used  in  this  thesis. 


Function  Subroutine  VALU 

The  VALU  function  is  called  by  the  DISSPLA  package  when 
a three  dimensional  plot  is  desired.  The  complex  aperture 
function  and  intensity  data  are  always  stored  in  the  complex 
array  DATA,  however,  DISSPLA  plots  only  real  data.  VALU 
assists  DISSPLA  by  providing  a real  array  of  data. 

SPECKLE  Main  Program 

This  section  covers  the  flow  of  data  through  the  SPECKLE 
program  as  outlined  in  Fig.  22.  The  operations  of  the  sub- 
routines and  the  input  data  formats  have  been  discussed  in 
previous  paragraphs  and  will  not  be  covered  in  detail  in  this 
section  of  the  thesis. 

The  first  operation  performed  by  SPECKLE  is  to  call 
subroutine  HEIGHS  which  will  input  the  normalized  surface 
height  measurements  and  the  other  program  variables.  Con- 
trol returns  to  the  main  program  where  the  normalized  sur- 
face height  data  is  output  on  the  line  printer  to  provide 
a numerical  record  of  the  input  data.  The  program  then 
checks  the  plotting  option  flags  (IN2D  and  IN3D)  and  deter- 
mines if  one  and/or  three  dimensional  plots  of  the  normal- 
ized surface  height  data  are  requested. 

SPECKLE  then  calls  COMAPE  to  calculate  the  complex 
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START 


PROGRAM  SPECKLE 


\ ^ 


CALL 

HPTnHc; 


LINE 

PRINTER 


User  supplied  subroutine 
Inputs ; 

I ROW,  ICOL,  MRDAT,  MCDAT , 
Plotter  requests 
Aperture  shape 
Normalized  height  data 


Output  normalized  height  data 


One  dimensional  plot  of 
height  data 


Three  dimensional  perspective 
plot  of  height  data 


I 


I 


Fig.  22  Flow  Chart  of  the  SPECKLE  Program 


(1  of  4) 
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r„,_. 


CALL 

COMAPE 


User  supplied  subroutine  to 
calculate  complex  aperture 
function 


CALL 

FOURT 


Call  FOURT  subroutine 


Calculate  intensity  field  from 
electric  field  distribution 


DATA  TO 
SECFFT 


Save  central  row  of  intensity 
pattern  for  frequency  content 
analysis 


NORMALIZE 

INTENSITY 


Normalize  intensity  field  i 

distribution  i 


PLOT 


One  dimensional  plot  of 
intensity  pattern 


Fig.  22  Flow  Chart  of  the  SPECKLE  Program  (2  of  4) 
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LINE 

PRINTER 


CALL 

FOURT 


NORMALIZE 

DATA 


LINE  i 

PRINTER 


Three  dimensional  perspective 
plot  of  intensity  distribution 


Output  central  row  of 
intensity  field  to  line 
printer 


Send  content  of  SECFFT  to 
FOURT  routine  for  frequency 
analysis 


Normalize  frequency  content 
data  for  plotting 


Output  normalized  central 
row  of  frequency  data 


Fig.  22  Flow  Chart  of  the  SPECKLE  Program  (3  of  4) 
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One  dimensional  plot  of 
frequency  data 


End  of  Program  SPECKLE 


Fig.  22  Flow  Chart  of  the  SPECKLE  Program  (4  of  4) 
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aperture  function.  The  resulting  data  is  passed  to  the 
FOURT  subroutine,  which  returns  the  two  dimensional  fast 
FOURIER  transform.  The  SHIFT  subroutine  then  completes  the 
symmetrical  rearrangement  on  the  data  to  obtain  the  zero 
frequency  magnit  ide  at  the  center  of  the  DATA  array.  The 
far  field  idffraction  pattern  is  finally  obtained  by 


DATAd , J) 


DATA(I,J) (Complex  Conjugate  DATA(I,J)) 
(2) (ALAMBA) ^ (RANGE) ^ 


(61) 


where  the  distance  between  samples  of  the  intensity  pattern 
in  the  observation  plane  is 


AMEOB  = IALAMBA) (RANGE)  (g2 

(IROW) (AMESC) 

This  completes  the  calculation  of  the  intensity  speckle 
(diffraction)  pattern,  in  the  far  field,  as  given  by  Eq.  47 
from  the  theory  chapter. 

The  maximum  value  of  the  speckle  pattern  is  located  and 
used  to  normalize  the  intensity  data  in  the  DATA  array. 

The  results  of  this  normalization  is  an  intensity  pattern 
with  a maximum  of  one.  Plotter  flags  (INT2D  and  INT3D) 
are  checked  to  see  if  one  dimensional  and/or  three  dimen- 
sional plots  of  the  speckle  pattern  are  desired.  In 
addition,  the  intensity  data  along  row  IROW/2  of  the  array 
DATA  is  stored  in  array  SECFFT  for  future  use  and  is  also 
printed  out  on  the  line  printer  for  a numerical  record  and 


69 


to  obtain  better  resolution  than  provided  on  the  plots. 

The  remaining  portion  of  the  program  computes  the  fre- 
quency content  and  proper  scaling  factors  for  the  center  row 
of  the  two  dimensional  speckle  pattern.  This  is  accomplish- 
ed by  first  sending  the  data  previously  saved  in  array 
SECFFT  to  FOURT  for  a one  dimensional  Fourier  transform  and 
then  performing  a single  row  shift  operation  (recalling 
that  SHIFT  is  used  on  two  dimensional  arrays) . The  pattern 
is  normalized,  plotted  if  requested,  and  a listing  of  the 
amplitudes  is  provided  on  the  line  printer.  The  program  is 
then  terminated. 

SPECKLE  Program  Deck  Structure 

Figure  23  shows  the  deck  structure  required  to  run  the 
SPECKLE  program  on  the  CDC  6600  computer  at  Wright-Patterson , 
AFB. 

SPECKLE  Program  Example 

For  this  example,  the  normalized  surface  height  data  was 
an  15x15  array  (square  aperture  shape)  and  the  DATA  array  was 
dimensioned  at  128x128.  The  laser  wavelength  was  6328  8, 
range  was  1000  meters,  and  the  distance  between  surface 
height  measurements  was  0.1  millimeters.  All  output  plots 
were  requested. 

Figures  24  through  28  are  the  five  plots  produced  by  the 
SPECKLE  program.  Figures  24  and  25  show  the  normalized  scat- 
tering surface  height  distribution  for  the  rough  surface, 
a one  dimensional  cut  and  a three  dimensional  perspective. 
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Normalize^.;  Surface  Height  Data 


Fig.  25  One  Dimensional  Cut  Through  Fig. 
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OBSERVATION  PLANE  MEASUREMENT  AT  4-94><10  METERS/UNIT 


3-D  NORMRLIZED  INTENSITY  DISTRIBUTION 
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Fiq.  27  Throe  Diraeesional  rerspcct-ive  of  Speckle  rattern 


Figures  26  and  27  show  the  speckle  (diffraction)  pattern  at 
the  observation  plane.  Figure  26  shows  a one  dimensional 
cut  through  the  center  of  the  intensity  pattern,  plotted  on 
a log-linear  graph.  Figure  27  is  a three  dimensional  per- 
spective of  the  intensity  pattern,  plotted  with  three  linear 
axis.  Figure  28  is  a one  dimensional  plot  of  the  Fourier 
transform  of  Fig.  26.  By  noting  on  Fig.  28  that  the  normal- 
ized frequency  content  goes  to  zero  at  15.0  units  or  2.37 
cycles/meter  (15.0  units  x .158  scale  factor),  the  typical 
width  of  the  intensity  lobes  in  the  intensity  pattern  should 
be  0.422  meters  (1/frequency).  On  Fig.  26,  0.422  meters  is 
equivalent  to  8.54  units  ( 0 . 422/4 . 94E-2 ) , which  is  the  width 
of  a typical  speckle  lobe. 

Recalling  from  the  theory  chapter,  the  typical  diffrac- 
tion (speckle)  lobe  width  is  related  to  the  linear  size  of 
the  aperture  by 

Aperture  Width  = -^2 (63) 

Typical  Lobe  Width 

From  the  previous  paragraph,  the  typical  lobe  width  was 
0.422  meters.  Using  Eq.  63  the  aperture  width  should  be 
15E-3  meters.  This  was  the  size  of  the  reflecting  area. 

Conclusion 

The  SPECKLE  program  is  designed  to  calculate  and  plot 
the  speckle  pattern,  in  the  far  plane,  produced  by  rough 
surface  scattering  of  a laser  beam.  The  user  is  required  to 
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provide  two  subroutines,  one  which  inputs  the  normalized 
surface  height  measurements  to  the  program  and  one  which 
calculates  the  complex  aperture  function  model.  The  scat- 
tering aperture  cross  section  is  described  by  the  arrange- 
ment of  the  normalized  surface  height  measurements  in  the 
array  DATA  through  the  subroutine  HEIGHS.  The  profile  of 
the  incident  radiation  field  and  the  distribution  of  the 
surface  reflection  coefficients  can  be  modeled  in  the  com- 
plex aperture  function  subroutine  COMAPE.  One  model  for 
the  complex  aperture  function,  as  described  by  Goodman 
(REF. 6),  has  been  described  and  used  in  this  thesis,  however, 
COMAPE  enables  the  user  to  use  a different  model  if  desired. 

This  chapter  was  written  to  describe  in  detail  the 
operation  of  the  program  and  the  significance  of  each  var- 
iable. Appendix  A contains  a listing  of  the  SPECKLE  pro- 
gram and  the  two  user  supplied  subroutines  used  by  this 
author  to  analyze  speckle  patterns  produced  by  rough  surface 
(simulated  data)  scattering  of  a laser  beam. 


V.  Verification  of  the  Speckle  Program 

This  chapter  will  address  the  validity  of  the  computer 
model  used  to  generate  speckle  patterns  produced  by  rough 
surface  scattering  by  comparing  the  computer  simulations 
with  theory  and  experimental  results.  The  speckle  pattern 
phenomena  that  will  be  discussed  are  the  affects  on  typical 
speckle  size  from  range  and  scattering  surface  size  varia- 
tions and  the  changes  in  speckle  pattern  distributions 
caused  by  variations  in  scattering  surface  roughness  and 
surface  reflectivity.  Because  of  the  numerous  computer 
plots  required  to  completely  describe  the  various  character- 
istics, all  data  plots  appear  in  Appendix  B. 

Speckle  Patterns  from  Smooth  Surfaces 

The  simplest  diffraction  patterns  result  when  the  scat- 
tering surface  is  smooth,  i.e.,  no  phase  variations  across 
the  reflecting  area.  The  two  smooth  diffracting  apertures 
that  are  always  discussed  and  analyzed  in  optics  text  books 
(Refs. 7 and  12)  are  those  with  square  (rectangular)  and 
circular  cross  sections.  The  square  aperture  produces  the 
well  known  |sinc(x)  normalized  intensity  distribution  in 
the  far  field,  which  was  given  by  Eq.  48  and  repeated  here  as 


. ,irbx. 
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or 


I(x,y)  = 


sin  3' 

2 

'sin  Y 

6 J 

1 y 

(65) 


where  3 and  y represent  the  arguments  of  the  sine  functions 
in  Eq.  64.  Recalling  from  the  theory  section  that  X = laser 
wavelength,  Z = range,  and  b = aperture  width. 

The  extrema  of  I(x,y)  can  be  found  by  taking  the  deriv- 
ative of  Eq.  65  and  setting  it  equal  to  zero.  In  one  dimen- 
sion the  derivative  is 


dl(x)  _ 2 Sin  3 (3Cos3-Sin3) 


From  Eq.  66  the  irradiance  is  zero  (nulls)  when  3 = mr,  where 
n = 1,  2,  •••00.  The  maxima  are  found  from  the  solutions  of 
the  transcendental  function 

3 Cos  3 - Sin  3=0  (67) 

which  reduces  to 

Tan  3=3  (68) 

From  Eqs.  66  and  68  the  first  three  nulls  occur  at  3 = Tr , 2tt, 
and  3tt  and  the  first  three  maxima  occur  at  3 = 1.4303tt, 
2.4590TT,  and  3.4707tt. 

Figures  29  through  32  of  Appendix  B are  plots  produced 
by  the  speckle  model  program  for  a square  aperture,  simulated 
by  discrete  data,  where  A = 6328  8,  Z = 1000  m. , and 


80 


I 


i 


i 

1 

j 


i 


1 


b=1.5  millimeters.  Table  III  lists  the  values  of  the 
normalized  peaks  and  locations  of  the  nulls  as  predicted  by 


theory  and  calculated  by  the  program  as  shown  in  Fig.  31. 

I 

1 The  computer  listing  provided  when  Fig.  31  was  plotted  was 

i used  to  obtain  the  resolution  of  the  normalized  peaks  report- 


I 


ed  in  Table  III. 


It  should  be  noted  that  the  first  null  of 


; 


Fig.  31  does  not  go  to  zero.  The  reason  for  this  unexpected 
result  is  that  the  zero  value  did  not  fall  on  a discrete 
data  sample . 


Table  III 

Smooth  Surface  Square  Aperture 
Figure  31  (Appendix  B) 
Diffraction  Pattern  Characteristics 


Null 

Location 

from 

Theory  (m) 

Null 

Location 

from 

Program  (m) 

Peak 

Amplitude 

from 

Theory 

Peak 

Amplitude 

from 

Program 

First 

0.422 

0.420 

0.0472 

0.0448 

Second 

0.844 

0.840 

0.0164 

0.0167 

Third 

1.  266 

1,260 

0.0083 

0.0091 

Although  not  developed  in  the  theory  section,  the  nor- 
malized intensity  pattern,  in  the  far  field,  for  circular 
aperture  reflection  (Ref. 7: 352)  is 
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In  Eq.  69,  b is  the  radius  of  the  aperture,  x is  the  radius 
in  the  observation  plane,  is  a first  order  Bessel  func- 
tion of  the  first  kind,  and  3,  X,  and  Z are  defined  as 
before.  Equation  69,  when  plotted,  is  the  well  known  Airy 
function.  The  first  three  nulls  are  obtained  when  Jj  equals 
zero,  i.e.,  3 = 3.831,  7.016,  and  10.17  (Ref . 13 : 250) . From 
the  equation  of  3,  the  x values  can  be  calculated  if  X,  b, 
and  Z are  given.  The  derivative  of  Eq.  69  can  be  set  equal 
to  zero  to  obtain  the  maxima  values,  where  d (23) /d3  = 

J2  (23) • The  first  three  secondary  peaks  occur  at  3 = 5.135, 
8.417,  and  11.619  (Ref . 13 : 250) . 

Figures  33  through  36  of  Appendix  B are  plots  provided 
by  the  speckle  program  for  a simulated  discrete  circular 
aperture  where  X = 6328  8,  Z = 1000  m. , and  b = 1.5  millime- 
ters. Table  IV  lists  the  values  of  the  null  locations  and 
normalized  secondary  peak  amplitudes  for  the  data  shown  in 
Fig.  35  and  as  predicted  from  theory. 

Table  IV 

Smooth  Surface  Circular  Aperture 
Figure  35  (Appendix  B) 

Diffraction  Pattern  Characteristics 


Null 

Location 

from 

Theory  (m) 

Null 

Location 

from 

Program  (m) 

Peak 

Amplitude 

from 

Theory 

Peak 

Amplitude 

from 

Program 

First 

0.248 

0.247 

0.0175 

0.0140 

Second 

0.456 

0.454 

0.0042 

0.0037 

Third 

0.661 

0.667 

0.0016 

0.0012 
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Figures  29  through  36  and  Tables  III  and  IV  are  data, 
for  smooth  scattering  surfaces  (square  and  circular) , cal- 
culated from  theory  and  by  program  SPECKLE.  With  the  data 
already  analyzed  there  is  a good  correlation  between  theory 
prediction  and  computer  model  results.  The  small  differ- 
ences in  results  observed  in  Tables  III  and  IV,  can  be  con- 
tributed to  discrete  data  processing. 

Effect  of  Range  on  Typical  Speckle  Size 

The  typical  speckle  size  is  directly  proportional  to 
range  as  predicted  by  theory  and  as  shown  by  the  experimental 
investigation  of  speckle.  The  chapter  on  the  experimental 
investigation  of  speckle  patterns  shows  photographic  and 
photomultiplier  data  that  supported  the  theory. 

Figures  37  through  45  of  Appendix  B are  plots  provided 
by  the  SPECKLE  program  to  aid  in  the  analysis  of  the  effect 
of  range  variations  on  typical  speckle  size.  Figures  37  and 
38  are  plots  of  the  normalized  scattering  surface  height 
variations  and  Fig.  39  is  a three  dimensional  perspective 
plot  of  the  intensity  pattern  at  the  observation  plane. 
Although  these  plots  show  a rough  scattering  surface,  the 
effect  of  the  rough  surface  on  the  speckle  pattern  will  not 
be  discussed  in  this  subsection  of  the  chapter.  Figures  40 
through  45  show  one  dimensional  cuts  through  each  of  the 
three  dimensional  perspective  plots  of  the  intensity  pat- 
terns and  Fourier  transform  plots  of  the  respective  one 
dimensional  cuts.  For  these  calculations,  observation  planes 
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were  located  at  500  m. , 1000  m. , and  1500  m,  from  the  scat- 
tering surfaces  as  indicated  on  the  one  dimensional  intensity 
plots . 


On  Figs.  41,  43,  and  45  the  reference  (REF)  marks  indi- 
cate where  the  normalized  Fourier  transform  plot  goes  to 
zero.  The  inverse  of  each  reference  spatial  frequency  pro- 
vides a measure  of  the  typical  speckle  width  on  the  respec- 
tive one  dimensional  intensity  plot.  Table  V lists  the 
reference  spatial  frequencies  and  typical  speckle  sizes  for 
the  three  ranges  indicated.  From  the  data  in  Table  V it  is 
apparent  that  the  speckle  model  accurately  calculated  the 
effect  of  range  variations  on  typical  speckle  size. 


Table  V 

Typical  Speckle  Size  vs  Range 


Range  (m. ) 

500 

1000 

1500 

Reference  Frequency 
(Cycles/m. ) 

4.74 

2.37 

1.57 

Typical  Speckle  Size  (m.) 

0.211 

0.422 

0.633 

Plots  (Appendix  A) 

40  and  41 

42  and  43 

44  and  45 

Effect  of  Reflecting  Surface  Size  Variations  on  Typical 
Speckle  Size 

In  the  theory  section  it  was  shown  that  the  typical 
speckle  size  is  inversely  proportional  to  the  reflecting 
surface  size.  The  experimental  analysis  of  speckl-  r 
produced  by  rough  surface  scattering  supported  t 
through  photographic  and  photomultiplier  dat 
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From  Figs.  37  through  41,  and  46  through  50  of  Appendix  B, 
using  the  same  analysis  procedure  as  in  the  previous  section, 
it  is  shown  that  as  the  reflecting  surface  size  increases 
the  typical  speckle  size  is  reduced  proportionally.  Table 
VI  lists  the  results  from  the  analysis. 

Table  VI 

Typical  Speckle  Size  vs  Scattering  Surface  Size 

Figures  Surface  Linear  Typical  Speckle 

Measurement  (m.)  Size  (m.) 

37  through  41  15E-4  0.422 

46  through  50  31E-4  0.204 

Percent  increase  in  surface  size  +106.7 

Percent  decrease  in  typical  speckle  size  -106.8 

Equally  Spaced  Apertures  with  Zero  Reflectivity  Between 
Apertures 

Before  discussing  speckle  patterns  produced  by  scatter- 
ing surfaces  with  random  surface  height  variations,  pulsed 
aperture  shapes,  where  the  reflectivity  is  zero  between  aper- 
tures will  be  discussed.  The  term  "pulsed"  means  that  the 

aperture  is  a periodic  array  of  rectangular  or  circular  ’ 

reflecting  (transmitting)  surfaces.  The  areas  for  each  unique 
aperture  have  a reflection  (transmission)  coefficient  of 
100%  and  the  areas  between  apertures  a zero  reflection  (trans- 
mission) coefficient.  Figure  51  shows  a three  dimensional 

perspective  of  an  array  of  square  holes  illuminated  by  a j 

plane  wave  coherent  laser.  i 
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Representations  for  single  and  multiple  pulse  one  dimen- 
sional apertures  and  respective  intensity  patterns  are  shown 
on  Fig.  52  on  the  following  page.  The  single  pulse  discrete 
data  aperture  and  its  far  field  intensity  pattern  have  been 
discussed  throughout  this  thesis.  Recall  the  normalized  one 
dimensional  intensity  pattern  has  the  form  of 


I/Io  = 


^Sin(2^)l^ 


where  the  meanings  of  each  variable  are  given  on  Fig.  52. 


Fig.  51  Periodic  Array  of  Square  Apertures 
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Although  not  proven  here  (see  Ref. 7 or  Ref. 12),  the 
far  field  normalized  one  dimensional  intensity  pattern  pro- 
duced by  plane  wave  laser  illumination  on  a periodic  set  of 
square  apertures  is  given  by 


In  Eq.  71,  P is  the  number  of  apertures,  q is  the  width  of 
the  apertures,  b is  the  interval  between  apertures.  The 
bottom  drawing  of  Fig.  52  illustrates  a sequence  of  aper- 
tures where  P = 4.  Note,  c = q + b is  the  period  of  the 
sequence . 

Several  conclusions  can  be  made  by  analyzing  the  bottom 
drawing  of  Fig.  52  and  Eq.  71.  First,  there  are  P-2  minor 
lobes  with  nulls  positioned  at  multiples  of  Xz/Pq  between 
each  major  lobe  which  are  positioned  at  multiples  of  Xz/q. 
Second,  as  b becomes  smaller  and  q increases,  the  secondary 
major  lobe  amplitudes,  positioned  at  multiples  of  Xz/q, 
decrease.  And  finally,  in  the  limit  as  q increases  and  b 
goes  to  zero,  Eq.  71  reduces  to  Eq.  70. 

A sequence  of  reflecting  apertures,  with  zero  reflec- 
tion between  apertures,  Fig.  53  of  Appendix  B,  was  formu- 
lated to  test  the  effectiveness  of  the  speckle  program  to 
calculate  the  far  field  intensity  pattern  as  predicted  by 
Eq.  71.  Note,  the  normalized  amplitudes  of  the  pulses  were 
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intentionally  to  keep  the  exponential  term  in  the  complex 


aperture  model  (Eq.  59)  from  going  to  1.  Thus,  the  complex 
aperture  model  returned  a complex  set  of  data  to  the  SPECKLE 
program.  The  plots  returned  from  the  SPECKLE  program  are 
shown  in  Figs.  53  and  56  of  Appendix  B,  with  the  following 
analysis  results: 


a. )  Figure  54  shows  there  are  five  pulses  sepa- 

rated by  zero  reflection  areas.  From  theory 
there  should  be  three  minor  lobes  between 
major  lobes  in  the  one  dimensional  intensity 
pattern.  Figure  56  shows  the  correct  results. 

b. )  The  nulls  in  the  intensity  pattern  associated 

with  the  pulse  separation  of  0.3  mm.  (Fig.  54) 
should  be  at  ±42.7  units  on  Fig.  56.  The 
correct  results  are  obvious  on  Fig.  56. 

c. )  The  first  nulls  in  the  intensity  pattern  asso- 

ciated with  a total  aperture  width  of  3.2  mm., 
five  pulses,  (Fig.  54)  should  be  at  ±4  units 
on  Fig.  56.  Again,  the  program  calculated 
the  nulls  correctly. 


For  clarity,  the  |sinc(x)|^  curves  on  the  intensity  pattern 
associated  with  the  total  aperture  width  and  the  pulse  sepa- 
ration interval  are  indicated  on  Fig.  56.  The  program  cal- 
culated and  plotted  the  DFT,  as  predicted  by  theory,  the 
diffraction  pattern  for  a series  of  pulse  shaped  apertures 
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and  zero  reflection  areas  as  depicted  in  Fig.  51. 

Reflecting  Surfaces  with  Arbitrary  Height  Variations 

In  the  previous  example,  a scattering  surface  was  repre- 
sented by  a series  of  pulses  of  equal  width,  spacing,  ampli- 
tude, and  areas  of  zero  reflectivity  (Fig.  53) . That  data 
was  used  to  show  the  effect  of  zero  reflection  areas  on  the 
speckle  pattern.  A rough  reflecting  surface  can  be  repre- 
sented by  varying  the  amplitudes  of  the  pulses  that  repre- 
sent the  scattering  surface  and  not  include  areas  of  zero 
reflection.  A three  dimensional  perspective  of  this  type  of 
reflecting  aperture  is  shown  in  Fig.  57  of  Appendix  B. 

Figure  58  shows  a one  dimensional  cut  through  Fig.  57,  how- 
ever, since  the  pulses  over  the  aperture  do  not  have  equal 
amplitudes.  Fig.  57  is  only  one  of  several  different  one 
dimensional  perspectives  across  the  aperture.  The  three 
dimensional  perspective  of  the  intensity  pattern  produced 
by  the  complex  aperture  is  shown  in  Fig.  5<^  and  a one  dimen- 
sional cut  through  the  intensity  pattern  is  shown  in  Fig.  60. 
Some  of  the  basic  characteristics  of  diffraction  theory 
previously  discussed  are  evident  in  the  two  dimensional 
intensity  plot  (Fig.  60) . First,  the  pulse  repetition  for 
this  aperture  (Fig.  57)  is  the  same  as  that  for  the  previous 
example  (except  there  are  no  areas  of  zero  reflection)  and 
the  intensity  diffraction  null  at  ±42.7  units  is  retained. 
Secondly,  the  diffraction  null  at  ±4  units,  associated  with 
the  total  aperture  width  of  3 .2  millimeters  is  retained. 
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Finally,  there  are  five  distinct  pulses  across  the  aperture, 
one  representation  is  shown  in  Fig.  58,  and  again  as  in  the 
previous  set  of  data  there  are  three  minor  lobes  positioned 
between  each  secondary  major  lobe  in  the  diffraction  pattern 
(Fig.  60) . The  major  change  in  this  diffraction  pattern 
(Fig.  60) , compared  to  the  diffraction  pattern  of  the  pre- 
vious example  (Fig.  56) , is  the  reduced  amplitudes  of  the 
secondary  major  and  minor  lobes. 

To  extend  the  simulated  aperture  representation  closer 
to  that  of  a rough  surface,  three  additional  parameter  var- 
iations were  added  to  the  pulsed  scattering  surface  of  the 
last  example;  random  pulse  width,  random  pulse  repetition 
interval,  and  sloped  pulse  edges.  A three  dimensional  per- 
spective of  this  type  of  aperture  is  shown  in  Fig.  61  of 
Aopendix  B.  Figure  62  shows  a one  dimensional  cut  through 
Fig.  61,  however,  as  stated  before  this  is  only  one  of  many 
cuts  where  each  would  reveal  a different  one  dimensional 
perspective.  Figure  63  shows  the  three  dimensional  perspec- 
tive of  the  intensity  diffraction  pattern  and  Fig.  64  shows 
a one  dimensional  cut  through  Fig.  63.  The  most  obvious 
result  is  that  all  secondary  peak  amplitudes  are  greatly 
reduced,  when  compared  to  the  two  previous  examples.  Other 
phenomena  are : 

a.)  There  are  three  minor  lobes  positioned  between 
the  central  lobe  and  the  first  major  secondary 
lobe  positioned  at  +18  units  on  Fig.  64.  This 
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was  expected  since  there  are  five  pulses 
across  the  aperture,  although  they  are  not  of 
equal  width  or  separation, 
b.)  Figure  64  shows  one  minor  lobe  positioned 

between  two  major  secondary  lobes  at  ±18  and 
±27  units.  A first  intuition  would  be  to 
assume  that  some  portion  of  the  three  dimen- 
sional complex  aperture  was  being  represented 
as  three  pulses,  which  would  produce  one  minor 
lobe  between  two  major  lobes.  However,  since 
the  intensity  pattern  is  the  two  dimensional 
Fourier  transform  of  the  complex  aperture  func- 
tion this  simplified  generalization  cannot  be 
assumed. 

The  next  example  is  a scattering  surface  with  random 
height  variations  as  shown  in  Fig.  65  of  Appendix  B.  Figure 
66  shows  a typical  one  dimensional  cut  through  Fig.  65. 
Obviously,  each  unique  cut  through  Fig.  65  will  produce  a 
different  one  dimensional  plot.  Figure  67  is  the  three  dimen 
sional  persepctive  of  the  intensity  pattern  produced  by  the 
complex  aperture  function  model  and  Fig.  60  is  a one  dimen- 
sional cut  through  Fig.  67.  Figure  68  shows  that  there  is 
a random  distribution  of  secondary  major  and  minor  peaks.  A 
Fourier  transform  plot.  Fig.  69,  of  the  one  dimensional  inten 
sity  plot.  Fig.  68,  reveals  the  spatial  frequencies  present 
in  the  diffraction  pattern.  The  frequency  reference  marks 


(REF)  on  Fig.  68  are  related  to  the  reference  (REF)  marks 
on  Fig.  69  through  Table  VII  and  Eq.  72. 

Intensity  plot  units  = (Frequency  plot  units  x 
Frequency  plot  scale  factor  x Intensity  plot 

scale  factor)”^  (72) 

Since  31  units  on  the  frequency  plot,  Fig.  69,  is  equal  to 
4.9  cycles/meter,  this  converts  to  0.2  m.  for  a typical 
speckle  width.  Fig.  68.  Transforming  back  to  the  scattering 
plane,  this  data  gives  an  aperture  width  of  3.1  milimeters. 
This  is  the  correct  result  as  predicted  by  theory. 

Table  VII 

Diffraction  Pattern  Frequency  Content 

Figure  69  Figure  68 


REF  1 

(units) 

31 

4.1 

REF  2 

(units) 

18 

7 

REF  3 

(units) 

8.2 

15 

REF  4 

(units) 

3.7 

35 

REF  5 

(units) 

2 

64 

final 

example 

for  a circular 

rough  scattering 

is  shown  in  Figs.  70  through  74  of  Appendix  B.  Although  not 
discussed  in  this  thesis,  the  analysis  of  the  spatial  fre- 
quency content  of  the  intensity  plot  is  similar  to  that  for 
the  square  reflecting  surface  shape.  The  important  factor 
to  remember  when  analyzing  diffraction  patterns  produced  by 
circular  cross  sections  is  that  the  intensity  nulls  do  not 
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occur  at  integer  multiples  as  they  do  for  a diffraction 
pattern  produced  by  a square  scattering  surface. 

Summary 

The  SPECKLE  program  was  tested  by  first  simulating 
simple  apertures  that  produced  well  known  diffraction  pat- 
terns and  then  by  making  the  apertures  more  complex  to  repre- 
sent rough  surfaces.  The  program  calculated  and  displayed 
the  correct  far  field  diffraction  patterns  for  those  aper- 
ture shapes  with  well  known  results,  i.e.,  smooth  surfaces, 
square,  and  circular  reflecting  cross  sections.  Although 
the  intensity  patterns  could  not  be  compared  to  text  book 
results  for  rough  scattering  surfaces,  there  were  certain 
characteristics  that  were  checked  for  program  verification. 
These  characteristics  were  null  locations  associated  with 
the  total  aperture  width  and  a complex  distribution  of 
secondary  major  and  minor  lobes  associated  with  surface 
roughness  characteristics.  For  all  examples  the  program 
provided  correct  displays  and  accurate  numerical  results. 

It  was  shown,  through  example,  that  there  are  several 
ways  to  represent  the  profile  of  scattering  surface.  The 
electric  field  at  the  scattering  surface  after  reflection  is 
represented  by  the  complex  aperture  function  model 

E,(u,v)  = AIFIEL(u,v)REF(u,v)Exp(l^(u,v) 

« A 

l.+Cos (GLAZ) ) (73) 
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If  the  profile  of  the  incident  radiation  is  represented  by 
the  normalized  height  distributions  of  the  scattering  sur- 
face, h{u,v),  phase  variations  will  be  present  in  the  com- 
plex aperture  function,  Eg(u,v)  data.  When  pulses  are 
represented  in  the  reflecting  aperture,  a value  of  zero  for 
h(u,v)  will  not  represent  zero  reflection  between  pulses. 

As  seen  in  Eq.  73,  a zero  for  h(u,v)  forces  the  exponential 
to  one  and  E (u,v)  will  have  a real  value.  Areas  of  zero 

di 

reflection  must  be  represented  by  the  reflection  coefficient, 
REF(u,v),  for  correct  results. 

By  using  Eq.  73  it  is  possible  to  fill  the  entire  array 
DATA  with  normalized  surface  height  data  and  model  the  inci- 
dent radiation  cross  section  with  the  field  term,  AIFIEL(u,v) 
or  the  surface  reflection  coefficient  term,  REF(u,v).  By 
using  the  three  major  variables  in  Eq.  73,  AIFIEL(u,v), 
REF(u,v),  and  h(u,v),  it  is  possible  to  use  the  same  set  of 
normalized  su'-face  height  data  to  model  complex  reflecting 
scenarios.  As  an  example,  consider  a reflecting  surface 
with  normalized  height  data,  h(u,v),  that  is  randomly  pointed 
with  an  absorption  material  resulting  in  a complex  distri- 
bution for  REF(u,v).  The  incident  radiation  distribution, 
AIFIEL(u,v),  can  now  be  varied  to  represent  different  shaped 
laser  beams.  Thus,  the  speckle  model  provides  numerous 
options  for  investigating  far  field  speckle  patterns. 
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VI.  Thesis  Conclusion  Applicable  to  ALOT  System  Design 

This  final  chapter  is  provided  for  the  ALOT  systems 
engineers,  at  the  AFWL,  who  initiated  this  thesis  project. 

The  major  areas  that  were  to  be  investigated  were  typical 
speckle  size  in  a far  field  diffraction  pattern  and  induced 
spatial  frequencies  caused  by  a speckle  pattern.  These  char- 
acteristics for  speckle  patterns  have  been  answered  in  the 
thesis  and  are  summarized  in  the  following  paragraphs. 

The  typical  speckle  size  in  a far  field  diffraction 
pattern  is  determined  by  the  cross  section  of  the  laser  beam, 
range  from  target  to  detector,  and  laser  wavelength.  Equa- 
tion 74  gives  the  typical  speckle  size,  at  a far  field 
observation  plane,  produced  when  a plane  wave  coherent  laser 
beam  is  reflected  from  a square  rough  surface. 

Typical  speckle  size  = — — ■-r-  (74) 

^ Reflecting  width 

For  a circular  aperture  shape,  the  typical  speckle  size 
will  be  approximately  1.22  times  value  given  by  Eq.  74.  The 
1.22  factor  is  a characteristic  of  the  Airy  pattern  distri- 
bution produced  by  the  circular  aperture.  The  term  "approxi- 
mate" is  appropriate  because  the  intensity  lobe  widths  in 
an  Airy  pattern  are  not  equal  as  they  are  in  a |Sinc(x) 
pattern  distribution  produced  by  a square  aperture.  The 
highest  spatial  frequency  is  the  inverse  of  the  typical 
speckle  size. 
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It  should  be  noted  that  variations  in  reflectivity  and 
surface  roughness  introduce  secondary  frequencies  in  the  dif 
fraction  pattern  which  are  inversely  proportional  to  the 
characteristic  length  of  the  variation.  That  is,  the  high- 
est spatial  frequency  is  inversely  proportional  to  the 
total  reflecting  surface  size  but  lower  spatial  frequencies 
result  from  unique  characteristic  lengths  of  the  reflecting 
area . 
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Appendix  A 
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SPECKLE  Program  Listing  (Page  2 of  40) 
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